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EARTHQUAKE  RESPONSE  OF  FIXED  OFFSHORE  STRUCTURES 

Owen  Major  Kirkley,  Ph.D. 

Department  of  Civil  Engineering 

University  of  Illinois  at  Urbana-Champaign,  1973 

The  response  of  fixed  offshore  structures,  such  as  those  used  for 
oil  drilling  and  production  purposes,  to  earthquake  excitation  is  investigated. 
An  approximate  method  of  computing  maximum  response  quantities  in  the  structure 
is  developed.  The  approximate  method  uses  superposition  of  classical  normal 
mode  responses  as  determined  by  response  spectrum  techniques.  The  structure 
is  modeled  as  a  lumped  parameter  system,  the  earthquake  is  represented  by  the 
deterministic  time  history  record  of  the  ground  acceleration,  and  the  hydro- 
dynamic  forces  are  considered  to  be  those  computed  from  Mori  son's  equation. 

A  single  degree  of  freedom  oscillator  with  nonlinear  drag  character- 
istics is  subjected  to  various  representations  of  ground  motion,  including 
earthquake  motions.  The  existence  of  a  high  frequency  tail  in  the  response 
spectrum  is  shown,  and  a  method  for  linearizing  the  nonlinear  hydrodynamic 
drag  force  is  developed  and  shown  to  be  satisfactory.  A  method  is  developed 
for  conservatively  estimating  the  magnitude  of  equivalent  linear  hydraulic 
damping  in  the  system  based  on  velocity  characteristics  of  the  ground  motion. 
Rules  are  developed  for  the  construction  of  idealized  response  spectra  for 
several  different  response  quantities.  The  effects  of  the  nonlinear  hydro- 
dynamic  drag  force  on  the  inelastic  response  of  the  single  degree  of  freedom 
oscillator  are  investigated. 

The  principles  and  approximate  methods  developed  during  the  study 
of  the  single  degree  of  freedom  system  are  applied  to  the  multi -degree  of 


freedom  model  of  the  actual  structure.  The  linearization  technique  is  ex- 
tended to  the  multi -degree  of  freedom  system  in  such  a  manner  as  to  insure 
the  existence  of  classical  normal  modes  in  the  equations  of  motion.  A  series 
of  multi -degree  of  freedom  systems  is  subjected  to  one  of  the  earthquake 
records  and  the  linearization  technique  is  shown  to  be  satisfactory  in  com- 
puting the  response  quantities  of  interest.  The  response  of  a  series  of 
offshore  towers  to  an  earthquake  record  is  computed  by  the  superposition  of 
modal  responses  as  determined  by  both  exact  spectra  and  idealized  spectra. 
The  approximate  method  considers  the  effects  of  the  structure-fluid  inter- 
action by  conservatively  estimating  the  amount  of  linear  hydraulic  damping 
developed  as  a  result  of  this  interaction. 


CHAPTER  1 
INTRODUCTION 
1.1  General 

The  design  of  fixed  offshore  structures  is  frequently  based  on  the 
static  forces  resulting  from  the  loads  imposed  on  the  structure  by  wind, 
waves,  ocean  currents,  and  earthquakes.  Design  on  the  basis  of  static  forces 
is  satisfactory  when  the  response  of  the  structure  is  not  magnified  by  the 
dynamic  nature  of  the  loads.  However,  the  characteristics  of  the  structure 
and  the  forces  acting  upon  it  are  often  such  that  the  response  of  the  struc- 
ture to  wind,  waves,  currents,  and  earthquakes  is  dynamic  and  magnified  greatly 
over  the  static  response.  In  these  cases,  design  should  properly  be  based  on 
the  dynamic  response  of  the  structure.  Various  methods  are  available  for  cal- 
culating the  dynamic  response  of  a  land-based  structure  to  earthquake  forces, 
but  similar  calculations  for  offshore  structures  are  complicated  by  the  pres- 
ence of  nonlinear  hydrodynamic  forces  acting  on  submerged  portions  of  the 
structure.  Further,  these  nonlinear  forces  are  dependent  not  only  on  the 
physical  characteristics  of  the  structure,  but  also  on  the  actual  response  of 
the  structure  to  the  forces  acting  on  it. 

Many  offshore  structures  are  currently  in  use,  with  the  most  notable 
use  being  that  of  oil  drilling  and  production  platforms.  While  platforms  in 
water  depths  of  up  to  300  feet  are  fairly  common  now,  the  search  for  new 
sources  of  oil  and  gas  may  result  in  many  offshore  structures  situated  in 
greater  water  depths.  The  location  of  oil  resources  off  the  coast  of  California, 
in  the  Gulf  of  Alaska,  and  other  portions  of  the  circumpacific  earthquake  belt 


as  well  as  other  regions  of  seismic  activity  emphasizes  the  need  to  be  able 
to  predict  the  response  of  offshore  structures  to  earthquake  motions.  Drill- 
ing and  production  platforms  represent  a  very  large  investment  in  cost,  and 
large  numbers  of  personnel  are  involved  in  the  operation  of  such  a  structure. 
Failure  of  one  of  these  platforms  could  result  in  many  lives  being  lost  as 
well  as  the  loss  of  investment  both  from  the  cost  of  the  platform  and  the  loss 
of  production  time. 

1.2  Description  of  the  Problem 

An  offshore  structure  of  the  type  mentioned  in  the  preceding  sec- 
tion is  shown  schematically  in  Fig.  1.1(a)  and  is  a  continuous  system  with 
an  infinite  number  of  degrees  of  freedom.  Such  structures  are  space  frames 
composed  of  circular  pipe  sections  generally  welded  together,  with  diagonal 
bracing  as  shown  in  the  figure.  The  top  portion  of  the  structure  may  consist 
of  more  than  one  level  and  contains  the  heavy  equipment  necessary  for  oil 
drilling  and  production  and  often  contains  facilities  for  berthing  and  messing 
of  the  crew.  The  operating  and  production  portion  of  the  structure  is  gener- 
ally quite  heavy.  The  research  reported  in  this  thesis  is  based  on  a  lumped 
parameter  model1  of  the  structure  with  mass,  stiffness,  forces  and  other 
characteristics  of  the  structure  lumped  into  a  discrete  number  of  nodes,  as 
shown  in  Fig.  1 .1 (b). 

The  equations  of  motion  for  the  system  shown  in  Fig.  1.1(b),  sub- 
jected to  ground  motion  in  quiescent  water,  are,  in  matrix  form: 

rmj{xs>  +  [Cs]{xs  -  xg}  +  [K]{xs  -  xg}  =  {P}      (1.1) 


* 

Numbers  refer  to  List  of  References  at  end  of  text. 


where 

Fm-J  =  mass  matrix  of  the  system 

[C  ]  =  structural  damping  matrix  of  the  system 

[K]  =  stiffness  matrix  of  the  system 

x  =  absolute  acceleration  of  node  j 

x      =  absolute  velocity  of  node  j 
Sj 

x      =  absolute  displacement  of  node  j 
Sj 

x      =  ground  velocity 
g 

x      =  ground  displacement 

The  symbol  [*  J  denotes  a  matrix  in  which  the  off-diagonal  elements  are  zero, 
No  structural  dampers  are  shown  on  the  model,  but  the  damping  matrix  [C  ]  can 
take  a  variety  of  forms  as  will  be  discussed  in  a  later  chapter. 

The  hydrodynamic  forces  acting  at  each  node  of  the  system  are  con- 
sidered to  be  those  computed  from  Morison's  equation2  for  circular  cylinders 
as: 


{P}  =  -  rmj{y  -  [M<DJ{xs|xs|}  (1.2) 


where 


["m  ~J  =  added  mass  matrix3  =  p["C,VJ 

p  =  mass  density  of  the  fluid  surrounding  the  structure 

C,  =  empirical  inertia  coefficient  at  node  j 
J 

V.  =  enclosed  volume  at  node  j 

["ICj  =  drag  factor  matrix  =  2"p["CJ\J 


Cn   =  empirical  drag  coefficient  at  node  j 
J 

A.   =  projected  area  at  node  j 

J 

If  the  coordinate  transformation,  {Z}  =  {x  -  x  },  is  made,  where 

^    y 

{Z}  is  a  vector  of  a  motion  quantity  of  a  node  relative  to  that  of  the  ground, 
Eqs.  1.1  and  1.2  can  be  combined  to  yield: 

[*M J{Z>  +  [CS]{Z}  +  [K]{Z>  =  -rMj{xg}  -  rKDJ{xs|xs|}  (1.3) 

where 

["Mj  =  virtual  (or  total)  mass  matrix  =  ["rnj  +  N  J 
x     =  ground  acceleration 

Equation  1.3  differs  from  the  equations  of  motion  for  a  land-based 
structure  subjected  to  earthquake  forces  in  two  important  aspects.  First, 
the  presence  of  the  added  mass  term  causes  the  structure  to  possess  natural 
frequencies  and  mode  shapes  which  are  different  from  those  of  a  structure  with 
the  same  stiffness  and  only  the  structural  mass  FmJ.  Secondly,  the  hydro- 
dynamic  force  term  ["KD^]{x  |x  1}  is  not  present  in  the  conventional  land-based 
earthquake  problem,  and  this  nonlinear  hydrodynamic  force  causes  the  major  dif- 
ficulty in  solving  Eq.  1.3. 

If  the  nonlinear  hydrodynamic  force  ["KD-J{x  |x  |}can  be  converted 
into  an  "equivalent"  linear  hydrodynamic  force,  say  ["-KnL-J{x  },  Eq.  1.3  can 
be  written  as 

^MJ{Z>  +  [Cs]{Z}  +  [K]{Z}  =  -FMj{xg}  -  rKDLJ{xs}      (1.4) 

Since  {Z}  =  {x  -  x  },  Eq.  1.4  can  be  rewritten  as 


FMJ{Z}  +  [Cs]{Z}  +  rKDL-J{Z}  +     [K]{Z} 

=     -^Mj{xg}  -   M<DLJ{xg}  (1.5) 

Equation  1.5  differs  from  the  equations  of  motions  for  a  land-based  multi- 
degree  of  freedom  damped  structure  of  mass  ["Mj,  subjected  to  a  base  exci- 
tation, in  the  linear  drag  force  D"  KDL-J{x  }. 

One  well-known  and  convenient  method  of  solving  simultaneous  linear 
differential  equations  of  motion  is  the  normal  mode  method.1  In  the  normal 
mode  method,  the  vector  of  a  response  quantity  expressed  in  one  coordinate 
system,  such  as  the  relative  displacement  between  each  node  and  the  base  {Z}, 
is  transformed  into  the  "normal"  coordinate  system  by  the  relation 

{Z}  =  [<f>]{Y}  (1.6) 

where 

[<J>]  =  a  matrix  of  the  natural  mode  shapes  of  the  system  with  the 

ith  column  in  the  matrix  representing  the  deflected  shape 

of  the  ith  mode  of  the  system, 
and      {Y}  =  a  vector  of  modal  response  quantities  with  the  ith  element 

in  the  vector  representing  the  modal  response  quantity  in 

the  ith  mode. 

Thus,  Eq.  1.6  expresses  the  response  at  each  node  in  the  model  as  a  combina- 
tion of  the  responses  at  that  node  in  each  of  the  natural  modes  of  vibration 
in  the  system.  The  natural  mode  shapes  in  the  modal  matrix  [(J)]  are  orthogonal1, 
with  respect  to  [*MJ  and  [K]  and  will  be  orthogonal  with  [C  ]  if  the  elements 
of  [C$]  are  specified  in  a  certain  way.21  Substitution  of  Eq.  1.6  into  Eq.  1.5 
yields, 


FMjfoHY}  +  [CS][4>]{Y>  +  r~KD|_J[(})]{Y}  +     [K][<J>]{Y} 

=     -rMj{xg}  -  rKDLJ{Ag}  (1.7) 

Pre-multi plying  Eq.   1.7  by  [<(>]     yields, 

[cf)]Tr-Mj[(J)]{Y}  +  [c})]T[Cs]U]{Y}  +     [4)]TrKDLJ[(t)]{Y} 

+     Dl>]T[K][<t>]{Y}  =     -[<f>]TrMj{xg}  -  U]TrKD|_J{xg}  (1.8) 

Recognizing  the  orthogonal   relationships  between  the  modal  matrix  [<}>]  and 
other  matrices  in  Eq.   1.8,  the  following  matrices  are  defined: 

*  T 

r*M  ^]  =  [cf>]  ["M -J[<})]         =     a  diagonal  modal  mass  matrix 

*  T 

[*C  J  =  M   [C  ][<(>]    =  a  diagonal  modal  structural  damping 

matrix.  The  form  of  [C  ]  which  permits 
this  is  discussed  later. 
C"K  J  =  [<f>]  [K][<J>]     =  a  diagonal  modal  stiffness  matrix 

[Km]  =  M   r*KnlJ[<j>]    =  a  modal  linear  hydraulic  damping  matrix 
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which  is  not  necessarily  diagonal 


r*w  J  =  ["M  -J"  [""K  -J  =  a  matrix  of  circular  natural  frequencies 

of  the  system 

Using  these  quantites,  Eq.  1.8  becomes 

C-M*.J{Y}  +  ^C*J{Y}  +  [K^_]{Y}  +  M(*J{Y} 

=     -[«D]T^Mj{xg}     -     U]TrKDLJ{xg}  (1.9) 


If  the  modal  linear  hydraulic  damping  matrix  [IC. ]  can  be  made  dia- 
gonal, then  the  set  of  simultaneous  differential  equations  of  motion  repre- 
sented by  Eq.  1.9  will  be  completely  uncoupled.  Suppose  that  the  nondiagonal 

matrix  [Km  ]  can,  in  some  satisfactory  manner,  be  converted  into  a  diagonal 

* 

matrix,  say  ["""Km  -J-  Then  Eq.  1.9  can  be  written 

FM*.J{Y}  +  rc*^m  +  ITKD*J{Y}  +  [M<*^{Y} 


=     -[>]TFMj{xg}     -     Df>]TrKDLJ{xg} 


(1.10) 


Or,  for  the  ith  mode, 


M..  Y.  +  C    Y.  +  Knl   Y.  +  K..  Y. 
n  l     s..  l      DL  .  i     n  i 
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where 


nn  =  number  of  nodes  in  the  model 


Equation  1.11  can  also  be  written  as 


Y,  +  23,  to..  Y.  +  2B.   to..  Y.  +  uj.  Y. 
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(1.12) 


where 


n 


■1   2H1fUl1 


=  structural  damping  ratio  in  the  ith  mode 


* 

KDL.. 
$.   =  — *   11    =  linear  hydraulic  damping  ratio  in  the  ith  mode 

•  •• 

If  the  modal  response  quantity  of  interest,  Y,  Y,  or  Y  can  be  ob- 
tained from  Eq.  1.12  for  each  mode,  then  the  total  response  of  the  structure 
can  be  obtained  by  using  the  relationship  given  in  Eq.  1.6.  The  modal  re- 
sponse quantities  vary  with  time,  and  the  time  variation  is  different  for 
each  mode,  so  an  exact  calculation  of  the  total  response  of  the  structure 
can  only  be  made  if  the  exact  time  histories  of  the  modal  responses  are  known 
and  superimposed  on  one  another.  Generally,  the  maximum  response  quantity 
at  each  node  of  the  model  is  desired,  but  the  modal  maxima  do  not  necessarily 
occur  at  the  same  time.  An  approximation  to  the  actual  maximum  response  at 
each  node  of  the  model  can  be  made  by  determining  the  maximum  response  in  each 
mode  and  then  combining  the  modal  responses  in  some  fashion.  The  response 
spectrum  is  a  convenient  way  to  determine  the  maximum  modal  response.11"13 

1.3  Summary  of  Previous  Research 

Equation  1.3,  or  a  similar  form,  has  been  solved  by  other  investi- 
gators using  various  methods.  Much  of  the  existing  work  in  this  area  has 
been  devoted  to  the  response  of  offshore  structures  to  wave  forces,  but  the 
nonlinear  difficulties  are  somewhat  similar  for  wave  forces  and  earthquake 
forces  acting  on  submerged  or  partially  submerged  structures.  Further,  the 
bulk  of  recent  investigations  has  been  devoted  to  the  use  of  stochastic  repre- 
sentations of  wave  and  earthquake  motions. 

Burke  and  Tighe1*  integrate  the  equations  of  motion  directly  for  a 
series  of  towers  subjected  to  different  deterministic  ground  acceleration 


records  by  using  an  iterative  technique  which  computes  the  response  quanti- 
ties of  interest  at  step-by-step  time  points  throughout  the  earthquake  record. 
For  a  tower  in  800  feet  of  water  subjected  to  four  accelerograms  (normalized 
to  the  intensity  of  the  1940  El  Centro  earthquake),  the  effect  of  the  water- 
structure  interaction  was  found  to  reduce  the  maximum  base  shear  and  the  maxi- 
mum overturning  moment  by  about  33  percent  from  the  values  found  when  the 
water-structure  interaction  was  not  included  in  the  solution.  This  method, 
while  permitting  the  direct  use  of  the  structural  properties  of  the  system, 
including  the  nonlinear  drag  force,  requires  the  use  of  a  digital  computer 
and  the  availability  of  actual  or  simulated  earthquake  records.  For  purposes 
of  preliminary  design  of  a  structure  for  a  "design"  earthquake  of  the  designer's 
choice,  however,  the  step-by-step  technique  is  costly,  time  consuming,  and 
requires  that  a  time  history  of  the  "design"  earthquake  be  available.  Pre- 
liminary design  would  best  be  served  by  a  rapid,  approximate  method  of  solu- 
tion which  permits  the  use  of  selected  parameters  for  the  design  earthquake 
and  accounts  for  the  effects  of  the  nonlinear  drag  forces,  yet  yields  satis- 
factory results. 

Penzien,  in  association  with  other  investigators,6"8  linearizes 
the  drag  forces  in  Eq.  1.3  by  replacing  the  nonlinear  force  vector  with  a 
linear  force  vector  which  performs  the  same  amount  of  work  during  the  history 
of  the  response.  This  linearization  results  in  an  "equivalent"  linear  hydraulic 
damping  matrix  acting  on  the  system  in  much  the  same  manner  as  the  structural 
damping  matrix  acts.  The  linear  equations  of  motion  are  then  transformed  into 
a  normal  coordinate  system  and  solved  by  the  normal  mode  method  of  superposition.1 
Unfortunately,  the  linear  hydraulic  damping  matrix  does  not  uncouple  when  the 
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equations  are  transformed  into  the  normal  coordinate  system.  As  a  result, 
an  uncoupling  (or  diagonal izing)  technique,  similar  to  the  original  linear- 
izing technique,  must  be  employed  in  order  to  separate  the  equations  into  the 
desired  modal  equations.  The  linearizing  and  diagonal izing  are  accomplished 
by  an  iterative  process  which  uses  actual  response  quantities  of  the  system 
as  they  are  developed,  and  the  procedure  is  time  consuming  and  requires  great 
care  in  its  formulation  and  use.  For  towers  subjected  to  a  stochastic  repre- 
sentation of  the  ground  motion,  equivalent  linear  hydraulic  damping  ratios  on 
the  order  of  9  percent  of  critical  are  reported  for  the  first  two  modes  of 
towers  in  deep  water,  with  "s/ery  much  smaller"  hydraulic  damping  ratios  for 
the  tower  in  shallow  water.6 

Foster9  solves  the  problem  of  an  offshore  tower  subjected  to 
stochastic  representation  of  wave  forces  by  linearizing  the  drag  force,  then, 
recognizing  that  nonproportional  hydraulic  damping  may  result,  solving  the 
linear  equations  of  motion  using  complex  eigenvalues  and  eigenvectors.  This 
procedure  recognizes  that,  unless  the  linearization  of  the  drag  force  is  done 
in  a  certain  way,  classical  normal  modes  may  not  necessarily  exist  in  the  equa- 
tions, thus  the  normal  mode  method  of  superposition  is  not  a  rigorously  cor- 
rect approach. 

Nath  and  Harleman10  conclude  that,  for  large  diameter  cylinders  in 
deep  water  waves,  the  nonlinear  drag  force  can  be  neglected,  and,  since  shallow 
water  structures  are  relatively  stiff,  and  therefore  the  dynamic  response  to 
waves  is  not  greatly  magnified  over  the  static  response,  design  of  shallow  water 
structures  can  proceed  on  a  static  basis.  These  conclusions,  however,  were 
arrived  at  by  considering  the  response  of  structures  to  waves,  which  generally 
have  much  longer  period  components  of  motion  than  do  earthquakes. 


11 


1.4  Objective  and  Scope 

The  objective  of  this  thesis  is  to  investigate  the  effects  of  the 
nonlinear  hydrodynamic  drag  force  on  the  earthquake  response  of  single  degree 
of  freedom  and  multi -degree  of  freedom  systems  and  to  develop  a  satisfactory 
method  of  linearizing  this  force  in  order  that  conventional  methods  of  solving 
linear  differential  equations  may  be  employed  in  the  solution  of  the  equations 
of  motion. 

A  further  objective  is  to  develop  an  approximate  method  for  the 
solution  of  Eq.  1.3  by  the  superposition  of  normal  mode  responses,  using 
response  spectra  to  estimate  the  modal  response  quantities  necessary.  A  rapid, 
yet  fairly  accurate  method  of  solution  is  desirable  in  the  preliminary  stages 
of  the  design  of  an  offshore  structure  in  order  that  the  designer  may  have 
an  estimate  of  the  response  of  a  tentative  design  to  one  or  more  design  earth- 
quakes of  his  specification.  It  should  be  kept  in  mind  that  the  approximate 
solution  developed  in  this  thesis  is  no  substitute  for  a  more  detailed  analysis 
of  the  structure  once  the  design  is  in  its  final  stages. 

In  Chapter  2,  the  effects  of  the  nonlinear  hydrodynamic  drag  force 
on  the  response  of  a  single  degree  of  freedom  oscillator  are  explored.  A 
reasonable  method  of  linearizing  the  drag  force  and  a  method  for  conservatively 
predicting  the  resultant  linear  drag  effects  are  developed  and  shown.  A  method 
is  developed  for  estimating  certain  response  quantities  of  the  oscillator  in 
order  to  determine  more  accurately  the  range  of  hydraulic  drag  and  inertia 
coefficients  which  should  be  used  in  the  force  equation.  The  response  of  the 
oscillator  in  the  inelastic  region,  assuming  elasto-plastic  force-displacement 
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characteristics  of  the  spring,  is  shown.  Finally,  the  construction  of  ideal- 
ized response  spectra,  which  include  the  effects  of  hydraulic  drag,  is  shown. 

The  concepts  developed  in  Chapter  2  for  the  single  degree  of  freedom 
system  are  extended  to  the  multi -degree  of  freedom  system  in  Chapter  3.  The 
nonlinear  drag  force  is  linearized  in  such  a  manner  as  to  insure  the  existence 
of  classical  normal  modes  in  the  equations  of  motion  in  order  that  normal  mode 
superposition  may  be  used  as  an  approximate  solution  to  the  problem.  The  method 
of  linearization  is  shown  to  be  accurate  for  a  series  of  five  degree  of  free- 
dom systems  which  are  representative  of  the  type  of  structure  being  investi- 
gated. A  means  of  estimating  the  amount  of  equivalent  linear  hydraulic  damping, 
using  ground  motion  velocity  characteristics,  is  developed  and  is  shown  to 
yield  good  results.  The  construction  of  an  idealized  response  spectrum  for 
each  mode  of  the  multi -degree  of  freedom  system  is  discussed  and  described, 
and  sample  problems  are  solved  by  the  superposition  of  normal  mode  responses. 

The  application  of  the  approximate  method  to  a  series  of  actual  off- 
shore towers  is  demonstrated,  and  the  results  of  calculations  of  the  response 
of  these  towers  are  discussed  in  Chapter  4. 

1 .5  Related  Topics 

In  this  section,  several  topics  of  general  interest  in  the  field  of 
earthquake  engineering  are  briefly  discussed  in  order  that  the  reader  may  more 
fully  understand  the  display  of  results  by  the  use  of  response  spectra,  the 
numerical  method  of  integrating  the  equations  of  motion,  and  the  characteris- 
tics of  the  earthquake  records  used  in  the  research. 
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1.5.1  Response  Spectra 

The  equation  of  motion  for  a  damped,  single  degree  of  freedom  os- 
cillator in  air,  subjected  to  a  base  disturbance  is 

mxs  +  Cs(xs  -  xg)  +  K(xs  -  xg)  =  0  (1.13) 

By  making  the  transformation  Z  =  x  -  x  ,  Eq.  1.13  becomes 

mZ  +  C$Z  +  KZ  =  -  mx  (1.14) 

Or,  in  an  alternate  form 

1     +  23$o)Z  +  w2  Z  =  -  x  (1.15) 

The  solution  of  Eq.  1.15  depends  on  the  characteristics  of  the  ground  motion 
function,  the  amount  of  damping  in  the  system  as  represented  by  the  dimension- 
less  parameter  $  ,  and  the  natural  frequency  of  the  system,  to.  The  frequency 
parameter  enters  into  the  solution  as  a  dimensionless  product  of  co  and  a 
characteristic  time  of  the  base  motion.13 

A  plot  of  the  maximum  value  of  a  response  quantity,  without  regard 
to  sign,  as  computed  from  Eqs.  1.13,  1.14  or  1.15,  as  a  function  of  the  fre- 
quency of  the  system  is  called  a  response  spectrum  for  that  quantity.  The 
response  spectrum  is  a  convenient  way  to  represent  various  maximum  response 
quantities  of  interest  as  a  function  of  the  variables  involved  in  the  solution. 
The  response  quantity  of  greatest  interest  is  usually  the  maximum  spring  dis- 
tortion, |Z|   ,  which  occurs  during  the  response  history  of  the  oscillator, 
although  other  response  quantities  are  sometimes  of  interest.13 

The  response  spectrum  is  conveniently  shown  on  tripartite  logarithmic 
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graph  paper  as  shown  in  Fig.  1.2.  The  grid  of  this  representation  is  con- 
structed such  that  the  abscissa  represents  the  frequency  of  the  oscillator, 
the  ordinate  represents  velocity  units  of  response,  the  diagonal  lines 
sloping  up  to  the  right  are  lines  of  constant  displacement  units,  and  the 
lines  sloping  down  to  the  right  are  lines  of  constant  acceleration  units. 
In  this  thesis,  unless  otherwise  indicated,  the  response  spectra  shown  are 
deformation  response  spectra  in  which  the  pseudo-velocity,  w|Z|    is  plotted 
against  the  undamped  frequency  of  the  system.  The  displacement  quantity 

shown  is  the  maximum  spring  deformation,  |Z|   ,  and  the  acceleration  quantity 

max 

?. 
is  the  pseudo-acceleration,  u>  |Z|   .  Response  spectra  for  a  damped  single 

degree  of  freedom  oscillator,  in  air,  are  shown  in  Fig.  1.2  for  the  El  Centro 
earthquake.  The  shape  of  the  spectra  is  typical  for  strong  motion  earthquakes 
of  the  type  used  in  this  thesis.  More  detailed  information  on  response  spectra 
can  be  found  in  the  literature.11"11* 

Since  the  response  spectra  for  strong  motion  earthquakes  exhibit  the 
characteristic  shape  shown  in  Fig.  1.2,  the  spectra  can  be  idealized  by  straight 
line  segments.  Amplification  factors,  which  relate  the  idealized  spectrum  to 
the  ground  motion  quantities  x     ,  x     ,  and  x     ,  have  been  developed 
for  a  range  of  damping  ratios.11*  The  ability  to  construct  an  idealized  repre- 
sentation of  the  response  spectrum  for  a  given  set  of  conditions  permits  an 
estimate  of  the  response  to  be  made  without  having  an  actual  earthquake  record 
available. 

1.5.2  Numerical  Integration  of  the  Equations  of  Motion 

The  equations  of  motion  for  a  system  subjected  to  a  base  disturbance 


15 


can  be  solved  by  numerical  integration.  In  this  thesis,  numerical  integra- 
tion of  the  equations  of  motion  is  performed  by  the  Newmark  "Beta  Method."5 
In  the  calculations,  3  is  set  to  be  equal  to  1/6,  thus  representing  the  ground 
response  acceleration  as  linear  between  adjacent  time  points.  The  time  inter- 
val used  is  1/20  of  the  natural  period  for  frequencies  up  to  10  cps  and  1/40 
of  the  natural  period  for  frequencies  equal  to,  or  greater  than  10  cps.  These 
time  increments  insure  good  accuracy,  stability  and  convergence  in  the  solu-' 
tion.  Convergence  of  the  solution  at  any  time  point  is  assumed  to  have  oc- 
curred if  two  successive  iterations  yield  values  of  acceleration  which  are 
within  0.1  percent.  In  the  case  of  multi-degree  of  freedom  equations  of 
motion,  the  time  increment  is  based  on  the  shortest  natural  period  in  the 
system,  and  iteration  is  carried  out  until  the  computed  acceleration  at  each 
node  is  within  0.1  percent  of  the  value  computed  on  the  previous  cycle.  In 
all  cases,  the  integration  of  the  equations  of  motion  is  carried  out  for  a 
period  of  time  equal  to  one-half  the  natural  period  (the  longest  natural 
period  in  the  case  of  multi-degree  of  freedom  systems)  after  the  base  distur- 
bance ends. 

1.5.3  Earthquake  Records  Considered 

Four  actual  earthquake  records  are  used  in  this  thesis.  The  records 
are  described  numerically  in  Table  1.1  and  the  time  history  of  the  ground  ac- 
celeration, velocity,  and  displacement  for  each  record  is  shown  in  Figs.  1.3 
through  1.6.  In  each  earthquake  record,  a  very   small  (less  than  .005g)  arti- 
ficial acceleration  point  has  been  generated  near  the  end  of  the  record  so 
that  the  velocity  of  the  base  motion  is  zero  at  the  end  of  the  record. 
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Although  a  parabolic  baseline  correction16  had  previously  been  applied  to  the 
acceleration  records,  a  slight  residual  velocity  still  remained  at  the  end  of 
the  record  and  it  is  this  residual  velocity  which  is  removed  by  the  artificial 
acceleration  point.  The  four  records  chosen  are  considered  to  be  representa- 
tive of  strong  motion  earthquakes  on  firm  soil. 

1.6  Notation 

As  much  as  possible,  the  same  symbols  are  used  for  both  the  single 
degree  of  freedom  oscillator  and  the  multi -degree  of  freedom  system,  with 
matrix  and  vector  notation  indicating  the  multi-degree  of  freedom  system. 
Symbols  are  defined  as  they  first  appear  in  the  text.  For  ease  of  reference 
the  symbols  used  are  summarized  in  this  section. 

["A J  projected  area  matrix 

Cp  empirical  dimensionless  drag  coefficient 

Cj  empirical  dimensionless  inertia  coefficient 

[C  ]  structural  damping  coefficient  matrix 

L  C  J  modal  damping  coefficient  matrix  =  [<j>]  [C  ][<j>] 

^  s 

D  diameter  of  a  single  degree  of  freedom  oscillator 

{HDPF}  vector  of  participation  factors  for  hydrodynamic 

parameter  response  quantities 
{HDPF}  =  rK*^]"1[(f,]T[K]{I} 

{1}  unit  vector 
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[K]  structure  stiffness  matrix 

r~K  -J  modal  stiffness  matrix  =  [<{>]  [K][(j>] 

[^KnJ  drag  factor  matrix  =  2-p["CJ\J 


{Kn}  modal  drag  factor  vector  for  very  high  frequency 

system  =  ["M*^!'1^]1^  KnJ{I} 


VD' 


<\j 


* 


T)" 


["KL.-J  linear  drag  factor  matrix 


[Kn,  ]  modal  linear  drag  factor  matrix  =  [<j>]  [""IC. -J[<}>] 

9  max  dimensionless  drag  parameter 

["Mj  total  structure  dynamic  mass  matrix  =  Fm-J  +  [~m  J 


J, 


C"M  J  modal  mass  matrix  =  I>]  FMJO] 

{P}  hydrodynamic  force  vector 

{PF}  participation  factor  vector  for  displacement  response 

C^VJ  enclosed  volume  matrix 

•       •• 

{Y},  {Y},  {Y}  displacement,  velocity  ,  and  acceleration  vectors 

{Qh  (Qh  {Q}  in  normal  coordinates 

{Z},  {Z},  {Z}  vectors  of  displacement,  velocity,  and  acceleration 

of  the  structure  relative  to  the  base 

Z  deformation  at  which  yielding  occurs 
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a  relative  acceleration  between  fluid  and  a  single 

degree  of  freedom  oscillator 

i  subscripts  indicating  mode  number 

j  subscripts  indicating  node  number 

["mj  structural  mass  matrix 

F"m  ^J  added  mass  matrix  =  pC^C.V^J 

nm  number  of  modes  in  the  multi -degree  of  freedom  system 

nn  number  of  nodes  in  the  lumped  parameter  model 

p(t)  hydrodynamic  force  acting  on  the  single  degree  of 

freedom  oscillator 

v  relative  velocity  between  fluid  and  a  single  degree 

of  freedom  oscillator 

{x  },  {x  },  {x  }       ground  displacement,  velocity  and  acceleration  vectors 
y    y    y 

xq  max'  xa  max'  xa  max   max"imum  displacement,  velocity,  and  acceleration  in 

the  ground  motion  record 

{x  },  {x  },  {x  }       vectors  of  structure  displacement,  velocity,  and 

acceleration 

y.j>  y.p  y.  displacement,  velocity,  and  acceleration  of  the  ith 

q.j>  q.j>  q.  mode  single  degree  of  freedom  equation  of  motion 

after  the  modal  participation  factor  has  been  taken 
out.  Expressed  in  normal  coordinates 
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a  constant  of  proportionality  between  the  linear  drag 

factor  ma- 


i 


factor  matrix  and  the  total  mass  matrix   ["Kpj.-J 


6.  linear  hydraulic  damping  ratio  in  the  ith  mode 


3  structural  damping  ratio  in  the  ith  mode 

1  C  * 

{y}  vector  of  modal  stiffness  factors  {y}   =  [(f)]  [K]{I} 


ductility  factor  for  inelastic  response 
'  'max 


zy 

M  mode  shape  matrix 

P  mass  density  of  the  fluid 


ok.  circular  natural  frequency  of  the  ith  mode 

o2j  =  rM  ^"Vk  ^ 
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CHAPTER  2 
SINGLE  DEGREE  OF  FREEDOM  OSCILLATOR 
2.1  General 

The  normal  mode  method  of  solving  Eq.  1.3  requires  the  solution 
of  single  degree  of  freedom  equations  of  motion.  For  this  reason,  and  to 
understand  more  fully  the  effects  of  the  nonlinear  hydrodynamic  drag  force 
on  the  response  of  a  simple  system,  Chapter  2  is  devoted  to  the  investigation 
of  the  effects  of  the  hydrodynamic  drag  force  on  the  response  of  a  single  de- 
gree of  freedom  oscillator. 

In  the  investigation  of  the  response  of  the  oscillator  with  non- 
linear drag  characteristics  it  is  shown  that  the  response  of  the  oscillator 
in  the  very  high  frequency  region  can  be  estimated  if  the  drag  properties  of 
the  system  are  known  and  if  the  maximum  velocity  of  the  base  motion  is  known. 
The  nonlinear  drag  force  is  linearized  by  equating  the  work  done  by  a  ficti- 
tious linear  force  with  the  work  done  by  the  real,  nonlinear  force.  This 
method  of  linearization  is  shown  to  be  accurate  until  the  high  frequency  re- 
gion of  the  spectrum  is  reached.  It  is  shown  that  the  linear  drag  force  is 
equivalent  to  a  linear  hydraulic  damper  in  which  the  force  is  dependent  upon 
the  velocity  of  the  mass,  and  that  the  magnitude  of  the  hydraulic  damping 
coefficient  is  dependent  upon  the  amount  of  structural  damping  in  the  system. 
A  means  of  estimating  the  amount  of  linear  hydraulic  damping  is  developed 
for  earthquake  motions,  and  the  estimate  is  shown  to  be  conservative  except 
in  the  low  frequency  range.  The  maximum  response  of  the  system  in  the  mid- 
frequency  range  to  earthquake  excitation  is  shown  to  be  dependent  upon  the 
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total  amount  of  linear  damping  present,  and  not  upon  the  ratio  in  which  the 
damping  is  distributed  between  structural  and  hydraulic. 

Response  quantities  which  may  be  used  in  order  to  select  reasonable 
values  of  the  inertia  and  drag  coefficients,  Cj  and  CD,  are  explored,  and  a 
method  of  estimating  these  response  quantities  is  developed.  The  inelastic 
response  of  an  oscillator  with  el  as to-plastic  force-displacement  character- 
istics of  the  spring  is  investigated  and  the  results  are  shown  in  the  form  of 
response  spectra. 

The  construction  of  idealized  elastic  response  spectra  which  include 
the  effects  of  the  nonlinear  hydrodynamic  forces  is  demonstrated  in  this 
chapter. 

2.2  Elastic  Deformation  Response  Spectra 

The  equation  of  motion  can  be  written  for  the  single  degree  of 
freedom  oscillator  as 


MZ  +  C$Z  +  KZ  =  -Mxg  -  KDxs|xs|  (2.1) 


or 


1    +  23  wZ  +  Jl     =  -x   -  -£- k   |x  I  (2.2) 

s  g     M   s1  s1  v   ' 


or 


Z  ♦  2B<i  +  U2Z  =  -«   -  (KD  *q  ""■)(  ag_)  J  |J 

y  g  max 


(2.3) 


As  pointed  out  in  Chapter  1,  the  solution  of  the  equation  of  motion  for  a 
land-based  single  degree  of  freedom  oscillator  (Eq.  2.1  without  the  nonlinear 
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hydrodynamic  drag  force  ICx  |x  |)  depends  on  the  characteristics  of  the  ground 
motion,  the  amount  of  damping  in  the  system  as  represented  by  the  dimensionless 
parameter  3  ,  and  the  natural  frequency  of  the  system,  co.  The  solution  of 
Eq.  2.1  depends  on  the  magnitude  of  KD  in  addition  to  these  parameters.  A 
dimensional  analysis  of  Eq.  2.1  yields  the  dimensionless  drag  parameter 

KD  xg  max 
K 

as  one  of  the  parameters  on  which  the  response  of  the  oscillator  depends. 

When  using  a  maximum  ground  acceleration  of  0.5  g,  values  of  0.01,  0.2,  and 

0.6  represent  a  reasonable  range  over  which  to  vary  the  dimensionless  drag 

parameter.  For  the  maximum  ground  acceleration  in  the  El  Centro  earthquake 

record  (0.312  g),  values  of  0.006,  0.08,  and  0.50  for  the  parameter  yield 

slightly  greater  variations  in  the  magnitude  of  KD.  Obviously,  a  zero  value 

of  the  drag  parameter  represents  a  "no-drag"  system  and  increasing  values  of 

the  parameter  represent  increasing  magnitudes  of  KL  in  the  system  for  fixed 

values  of  K  and  k'  m  . 
g  max 

Equation  2.1  is  solved  numerically  for  the  steady-state  maximum 
response  of  the  oscillator  when  the  base  acceleration  is 

x   =  16.1  sin  2frt 

where 

2 
x   is  in  ft/ sec 

t   is  in  seconds 

Deformation  response  spectra  (in  which  <*>|Z|    is  plotted  against  frequency) 
are  plotted  in  Figs.  2.1  and  2.2  to  show  the  results  of  the  calculations. 
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In  order  to  investigate  the  response  of  the  oscillator  to  a  slightly 
more  complicated  representation  of  the  base  disturbance,  the  base  acceleration 
is  represented  as  a  six  term  cosine  series.  The  numerical  characteristics  of 
the  series  used  in  the  response  calculations  are  shown  in  Table  2.1.  The 
series  has  the  features  of  zero  acceleration,  velocity,  and  displacement  at 
the  beginning  and  end  of  the  base  motion.  Equation  2.1  is  solved  numerically 
for  the  maximum  response  of  the  oscillator  when  subjected  to  five  cycles  of 
this  base  disturbance,  and  the  results  are  shown  in  the  form  of  deformation 
response  spectra  in  Figs.  2.3  and  2.4. 

Equation  2.1  is  solved  numerically  for  the  maximum  response  of  the 
oscillator  when  subjected  to  the  El  Centro  earthquake  record.  The  results 
are  shown  in  the  form  of  deformation  response  spectra  in  Figs.  2.5  through 
2.7. 

A  study  of  the  spectra  presented  in  Figs.  2.1  through  2.7  reveals 

several  interesting  characteristics  which  are  best  discussed  in  terms  of 

frequency  ranges  of  the  oscillator. 

In  the  low  frequency  region,  the  maximum  spring  deformation,  |Z|   , 

max 

approaches  that  for  the  case  of  K  *x  max/K  =  0,  regardless  of  the  amount  of 
hydrodynamic  drag  in  the  system  for  the  range  of  drag  parameter  used.  This 
is  reasonable  to  expect,  since  the  magnitude  of  the  drag  force  is  a  function 
of  the  oscillator  velocity,  x  ,  and  very   soft  (low  frequency)  systems  exhibit 
little  velocity  response  of  the  mass  (x  ). 

In  the  mid-frequency  region,  where  spectra  for  the  conventional 
"no-drag"  oscillator  show  an  amplification  of  the  pseudo-velocity  over  the 
maximum  ground  velocity,  the  spectra  for  systems  with  drag  also  show  an  ampli- 
fication of  the  pseudo-velocity  over  the  maximum  ground  velocity.  However, 
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the  presence  of  the  nonlinear  hydrodynamic  drag  force  causes  a  reduction  in 
the  response  of  the  oscillator  from  the  response  of  the  "no-drag"  oscillator. 
The  amount  of  reduction  is  greater  for  greater  values  of  the  drag  parameter 
up  to  some  frequency.  This  reduction  would  suggest  that,  in  the  mid-frequency 
range  of  the  spectrum,  the  nonlinear  hydrodynamic  force  can  be  replaced  by 
an  equivalent  linear  hydrodynamic  force  which  will  effect  the  response  in 
much  the  same  manner  as  the  presence  of  structural  damping  effects  the  re- 
sponse. The  amount  of  reduction  in  response  caused  by  the  drag  force  is  also 
seen  from  the  spectra  to  be  dependent  upon  the  amount  of  structural  damping 
in  the  system.  This  is  easily  noted  by  comparing  Figs.  2.5  and  2.6.  In  Fig. 
2.5,  the  system  has  no  structural  damping,  and  a  slight  amount  of  drag  (repre- 
sented by  ICx    /K  =  0.006)  is  seen  to  reduce  significantly  the  maximum 
J     D  g  max        '  a       J 

response  of  the  oscillator  in  the  frequency  range  of  about  1  to  10  cps.  In 
Fig  2.6  the  system  has  2  percent  structural  damping,  and  the  drag  parameter 
of  0.006  is  seen  to  cause  very   little  reduction  in  the  response  spectrum  be- 
tween 1  and  10  cps.  The  effects  of  structural  damping  on  the  amount  of  re- 
duction in  response  caused  by  the  hydrodynamic  drag  are  to  be  expected,  since 
the  drag  force  is  a  function  of  the  oscillator  velocity,  x  ,  and  this  velocity 
is  a  function  of  the  amount  of  structural  damping  in  the  system. 

In  the  high  frequency  region,  response  spectra  for  "no-drag"  oscil- 
lators typically  show  the  pseudo-acceleration  approaching  the  maximum  ground 
acceleration,  x     .  Response  spectra  in  the  high  frequency  region  for  the 
oscillators  with  nonlinear  hydrodynamic  force  exhibit  entirely  different 
characteristics,  however.  In  all  the  response  spectra  shown,  a  frequency  is 

reached  at  which  the  maximum  response  quantity  |Z|    for  the  oscillator  with 

1  'max 
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drag  exceeds  that  for  the  "no-drag"  oscillator.  In  fact,  for  the  systems 
with  larger  values  of  the  drag  parameter,  the  quantity  |Z|    is  seen  to  ap- 

1 1  id  a 

proach  a  relatively  constant  value.  The  magnitude  of  this  high  frequency 
"tail"  in  the  spectrum  can  be  estimated  as  follows: 

Consider  an  undamped  system  with  spring  stiffness  so  great  that 


s     g  s     g 


The  forces  on  the  oscillator  are 


F   .    =  F.   ..   +  F. 
spring     inertia     drag 


F         Mx   +  Knx  |x  I 
spring  =    s     D  s1  s1 


or, 


Fspring  "  %  *  KD(*/ 
Since  the  drag  parameter 

Kn  x 
D  q  max 

K 

is  held  constant,  KQ  must  increase  as  the  spring  stiffness  K  increases,  so 
a  frequency  will  be  reached  at  which  the  inertia  force  can  be  neglected  with 
respect  to  the  drag  force.  If  the  inertia  force  is  neglected,  then 

(F     )     ~  K  (x     )^ 
v  spring 'max  ~   Dv  g  max' 

And 

(F  )  K  K     x  x 

|7|  =       spring  'max      _D  ,•  x2  ,  p     g  maxu  g  maxN         (9  Ax 

1    'max  k  ~    K    ^  9  max;  l         K         Mx  ;         U*4; 

1X1X3  *  g  max 
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A  comparison  of  |Z|    obtained  at  high  frequencies  by  the  numerical 
integration  of  Eq.  2.1  with  |Z|    computed  from  Eq.  2.4  is  shown  in  Table 
2.2  for  all  systems  and  base  disturbances  considered.  A  study  of  the  data 
presented  in  Table  2.2  shows  that,  particularly  for  the  earthquake  records, 
the  high  frequency  tail  of  the  deformation  response  spectrum  can  be  satisfac- 
torily estimated  by  the  use  of  Eq.  2.4. 

When  an  acceleration  record  is  integrated  numerically  in  order  to 

obtain  the  velocity  and  displacement  records,  the  integration  is  done  from 

time  point  to  time  point  in  the  record.  This  double  integration  yields  the 

values  of  x   ,  =  12.79  in. /sec  and  x ,„  =  8.24  inches  for  the  El  Centro 

g  max  g  max 

record  as  shown  in  Table  1.1.  However,  if  the  acceleration  record  is  inte- 
grated to  find  the  ground  velocity  and  displacement  at  some  time  point  which 
is  between  time  points  in  the  acceleration  record,  linear  interpolation  of 
the  acceleration  record  must  be  performed  in  order  to  obtain  the  acceleration 
quantity  necessary  to  compute  the  ground  velocity  and  displacement  at  the 
intermediate  time  point.  This  interpolation  introduces  small  errors  into 
the  acceleration,  which  in  turn  introduces  slightly  larger  errors  into  the 
ground  velocity  which  in  turn  introduces  even  larger  errors  into  the  ground 
displacement  computed  at  the  intermediate  time  point.16  It  was  found  that, 
when  using  the  very  small  time  increments  necessary  for  stability  and  con- 
vergence in  the  numerical  integration  of  the  equation  of  motion  at  high  fre- 
quencies (above  about  8  cps),  the  interpolation  of  the  El  Centro  acceleration 

record  introduced  errors  into  the  magnitude  of  x  such  that  x     was  found 

g  g  max 

to  be  about  14.0  in. /sec  in  lieu  of  the  value  12.79  in. /sec  previously  ob- 
tained, even  though  double  precision  was  used  in  all  calculations.  In  order 
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to  achieve  comparability  of  results  in  the  comparison  of  |ZL_„  estimated  from 

max 

Eq.  2.4  and  |Z|    obtained  at  high  frequencies  from  the  integration  of  the 

Hlu  A 

equation  of  motion,  x   a  =14.0  in. /sec  is  used  in  Eq.  2.4  to  obtain  the 

g  max 

results  shown  in  Table  2.2  for  Knx„  m  /K  =  0.006  and  0.08. 

u  g  max 

2.3  Linearization  of  the  Drag  Force 

It  was  seen  in  Chapter  1  that  the  use  of  the  normal  mode  method  to 

solve  Eq.  1.3  depends  on  the  ability  to  linearize  the  nonlinear  hydrodynamic 

force  vector  r"Kn-J{x  |x  1}  in  a  satisfactory  manner.  The  linearization  of 
u    s '  s ' 

the  hydrodynamic  drag  force  in  the  single  degree  of  freedom  system  is  inves- 
tigated in  this  section.  It  is  desirable  to  replace  the  nonlinear  drag 
force,  ICx  |x  |,  in  Eq.  2.1  by  a  linear  drag  force,  say,  IC.  x  ,  which  will 
yield  the  same  maximum  response  in  the  system.  Consider  the  work  done  by  the 
two  forces: 15 


Work  done  by  the  nonlinear 

drag  force  over  a  period  of  time 


Work  done  by  a  linear  drag 
force  over  a  period  of  time 

Equating  the  amounts  of  work  done, 

•t 


t 


o 
t 


(Sxslxsl)(xs)dt 


(KDLxs)(xs)dt 


t    2 

KDxs|xs|dt 


KDLx2  dt 


KD<Asl^sl>       =       KDL<V 

•2,-    , 
/  <xs'x   ^    \ 
KDL     =     KD         ^2  (2*5) 

1     <xs>  ' 
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(The  symbol  <  >  denotes  the  average  value  of  the  quantity  enclosed  over  the 
time  of  the  response  of  the  oscillator.)  This  expression  is  the  same  expres- 
sion used  by  other  investigators6-9  to  linearize  a  nonlinear  force.  With  the 
nonlinear  drag  force  replaced  by  an  "equivalent"  linear  drag  force,  Eq.  2.1 
becomes, 

MZ  +  CsZ  +  KZ  =  -Mxg  -  KD|_xs 

MZ  +  C$Z  +  KD|_Z  +  KZ  =  -Mxg  -  KD[_xg  (2.6) 


or, 


where, 


Z  +  23  o)Z  +  23uo)Z  +  to2Z  =  -x   -  23,a)X  (2.7) 

s        h  g      h  g 


KDL 
B.   =  «j3-  =  equivalent  linear  hydraulic  damping  ratio 


It  now  remains  to  investigate  how  the  maximum  response  in  the  equi- 
valent linear  system  compares  with  the  maximum  response  in  the  actual  nonlinear 
system. 

Equation  2.6  is  solved  numerically  for  |Zl  „  using  a  linear  drag 

max 

factor  IC.  which  is  linearized  in  accordance  with  Eq.  2.5  for  the  single  de- 
gree of  freedom  oscillator  subjected  to  the  series  representation  of  the  base 
disturbance.  The  quantities  <x  | x  | >  and  <x  >  necessary  for  the  linearization 
were  obtained  from  the  integration  of  the  nonlinear  equation  of  motion  for  the 
same  set  of  conditions.  The  results  are  plotted  on  Figs.  2.8  and  2.9,  and  it 
is  seen  that,  for  this  base  disturbance,  the  linearization  is  \/ery   good  until 
the  high  frequency  tail  begins  to  appear.  This  might  have  been  expected, 
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since  a  linearization  of  a  weak  nonlinearity  would  normally  be  better  than  a 

linearization  of  a  strong  nonlinearity,  and  the  high  frequency  tail  is  the 

result  of  the  nonlinear  hydrodynamic  drag  force  becoming  quite  strong  and 

dominant  in  Eq.  2.1 . 

A  similar  set  of  calculations  is  made  for  the  El  Centro  earthquake, 

The  response  spectra  for  IC  and  IC.  are  shown  in  Fig.  2.10  for  3  =0.02  and 

KJ<   ,  /K  =  0.50.  It  is  seen  from  this  figure  that  the  linearization  for 
D  g  max  3 

the  earthquake  motion,  while  not  quite  as  good  as  the  linearization  for  the 
series  representation  of  the  base  motion,  is  still  reasonably  good  up  until 
the  high  frequency  tail  of  the  spectrum  is  reached.  A  comparison  of  |Z| 
computed  for  the  nonlinear  drag  system  with  |Z|    for  the  linear  drag  sys- 
tern  is  shown  in  Figs.  2.11  through  2.13  for  all  of  the  variables  considered. 
It  is  seen  from  a  study  of  these  figures  that  the  maximum  response  of  the 
linear  system  is  close  to  the  maximum  response  of  the  nonlinear  system  until 
the  high  frequency  region  is  reached.  The  figures  further  indicate  that  the 
linearized  response  is  slightly  on  the  conservative  side  of  the  actual  re- 
sponse until  the  high  frequency  region  is  reached.  The  accuracy  of  the 
linearization  in  the  high  frequency  region  is  expected  to  become  poor  as  the 
nonlinearity  is  becoming  more  and  more  dominant  in  the  system,  and  thus  the 
linearization  becomes  less  and  less  equivalent.  However,  Eq.  2.4  permits  a 
fairly  good  estimate  to  be  made  of  the  maximum  relative  displacement  in  the 
high  frequency  region. 

Since  the  linearization  of  the  drag  force  in  Eq.  2.1  is  seen  to 
yield  a  slightly  conservative  maximum  response  at  frequencies  below  the  be- 
ginning of  the  high  frequency  tail  in  the  response  spectrum,  a  method  of 
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estimating  the  linear  drag  factor  IC,  from  the  known  nonlinear  drag  factor 
Kn  is  now  developed  for  earthquake  motions.  A  study  of  Eq.  2.5  shows  that 
two  velocity  quantities,  <x  |x  j>  and  <x  >  are  necessary  in  order  to  linearize 
the  drag  force  and,  since  these  quantities  are  actual  response  quantities 
of  the  nonlinear  system,  they  cannot  be  calculated  a  priori.  A  series  of 
computations  is  macfe  in  which  the  quantity 

•2i-  i 
<xslxsl> 

is  computed  for  the  actual,  nonlinear  oscillator  subjected  to  the  El  Centro 
earthquake.  The  results  are  shown  in  Figs.  2.14  through  2.16.  From  these 
figures  it  is  seen  that,  for  the  range  of  variables  studied,  if  the  quantity 

•2,-  , 

<xslxsl> 

<> 
is  estimated  as 

^•2i-  , 

<x  x  > 

g1  q1 

•2 

for  the  particular  earthquake  considered,  then  this  estimated  value  will  be 
less  than  the  actual 

^•2i«  , 

s1  s1 

for  all  frequencies  except  low  frequencies.  Using  the  value  of 
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•2,-  i 

<X«l\l> 

<x2> 

g 

in  Eq.  2.5  to  linearize  the  drag  term  results  in  a  value  of  IC.  which  is  less 
than  the  value  of  IC.  which  would  result  from  using  the  actual  oscillator 
velocity  quantities.  When  this  lower  value  of  IC.  is  converted  into  the 
linear  hydraulic  damping  ratio,  8h>  shown  in  Eq.  2.7,  the  value  of  $.  is 
less  than  the  actual  value  of  $.  which  would  be  computed  had  the  actual  os- 
cillator velocity  terms  been  available  for  use  in  Eq.  2.5.  In  other  words, 
linearizing  the  drag  term  by  using 

•2i"  i 
<x  x  > 

g'g1 

-2 
<x  > 

g 

in  Eq.  2.5  leads  to  a  conservative  value  of  KD.  which,  in  turn,  leads  to  a 
conservative  value  of  the  maximum  response  in  all  except  the  low  frequency 
region  and  the  high  frequency  tail  region  of  the  spectrum.  The  unconserva- 
tive  nature  of  IC..  in  the  low  frequency  region  when  computed  from  the  ground 
velocity  terms  is  reasonable  to  expect,  however,  as  the  amount  of  hydraulic 
damping  is  dependent  on  the  velocity  response  of  the  oscillator,  and  wery 
soft  (low  frequency)  systems  have  little  velocity  response,  so  the  amount  of 
hydraulic  damping  should  approach  zero  as  the  system  becomes  softer  and  softer. 
Figures  2.14  through  2.16  further  show  that,  for  a  system  with  given  drag 
characteristics,  the  linear  drag  factor  KDL  estimated  by  using  the  quantity 

•2 1  •  i 

g'  g' 
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becomes  less  conservative  as  the  amount  of  structural  damping  increases  and 
as  the  amount  of  nonlinear  drag  increases.  The  linear  drag  factor  estimated 
in  this  manner  does  not  become  unconservative  (except  in  the  low  frequency 
region)  for  the  range  of  variables  studied. 

The  amount  of  conservatism  introduced  by  the  use  of 

•2,-  . 

g'  g 

in  linearizing  the  drag  factor,  IC,  is  seen  in  Figs.  2.17  through  2.19.  In 
these  figures,  the  spectrum  for  the  linear  drag  solution  is  seen  to  be  con- 
servative up  until  the  high  frequency  tail  begins.  The  linearized  spectra 
were  found  to  be  slightly  unconservative,  as  expected,  in  the  low  frequency 
region,  but  the  scale  of  the  figures  is  such  that  this  unconservatism  cannot 
easily  be  seen.  The  use  of  Eq.  2.4,  repeated  here  for  convenience, 

lZLax  ■  (T^  max>2  ■  <^H^>fe         (2.4) 
3  g  max 

is  also  seen  to  provide  a  yery   good  estimate  of  the  nonlinear  drag  spectrum 

in  the  high  frequency  region.  At  the  point  in  the  spectrum  at  which  the  high 

frequency  tail  just  begins  to  appear,  the  intersection  of  the  linear  drag 

spectrum  with  the  estimate  of  the  high  frequency  tail  provided  by  Eq.  2.4 

is  seen  to  be  slightly  unconservative.  The  amount  of  unconservatism  in  this 

transition  "saddle"  to  the  high  frequency  tail  is  reduced  as  the  amount  of 

nonlinear  drag  in  the  system  increases.  It  is  seen  from  the  figures  that, 

while  linearizing  the  drag  force  by  the  relation, 

•2,«  , 
<x  x  > 

KDL  =  KD  (    -29   }  <2-8> 

<xg> 
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yields  conservative  results  in  the  mid-frequency  region  of  the  spectrum,  the 

degree  of  conservatism  is  not  so  great  as  to  be  considered  too  conservative. 

The  presence  of  the  nonlinear  hydrodynamic  drag  force  has  been  seen 

to  cause  a  reduction  in  |Z|    from  that  computed  for  the  "no-drag"  oscillator 

1  'max 

at  a  given  frequency,  up  until  the  high  frequency  tail  region  is  reached.  A 
measure  of  how  much  of  this  reduction  is  accounted  for  in  the  computation  of 
| Z |    for  the  system  with  linear  drag  (linearized  by  the  ground  velocity) 
is  seen  in  Fig.  2.20.  The  ordinate  in  the  figure  is  the  ratio  of  the  reduc- 
tion (from  the  "no-drag"  case)  in  |Z|    when  the  actual  nonlinear  drag  force 

31  'max  3 

is  used  to  the  reduction  (again,  from  the  "no-drag"  case)  in  |Z|    when  the 

3  a/i  imax 

drag  force  has  been  linearized  by  the  ground  velocity  terms.  In  Fig.  2.20, 

values  of  the  ordinate  greater  than  1  represent  unconservative  results.  The 

unconservative  nature  of  the  linearization  in  the  low  frequency  region  and 

the  high  frequency  tail  region  is  noted.  While  the  linearization  of  low-drag 

systems  is  seen  to  account  for  less  than  50  percent  of  the  actual  reduction 

ln  ULav  ln  many  cases,  the  actual  reduction  for  low-drag  systems  is  quite 
max 

small  itself,  as  seen  in  Fig.  2.17.  It  is  further  seen  from  the  figure  that 
more  of  the  reduction  in  |Z|    is  accounted  for  in  the  high-drag  systems  than 
in  the  low-drag  systems. 

Figures  2.17  through  2.20  show  a  measure  of  the  amount  of  conserva- 
tism introduced  in  |Zl  v  by  using  the  ground  velocity  characteristics  to 
linearize  the  drag  force  when  the  structural  damping  ratio,  3  >  is  0.02. 
For  systems  with  less  structural  damping,  the  amount  of  conservatism  will 
be  slightly  greater,  and  for  systems  with  more  structural  damping  the  amount 
of  conservatism  will  be  slightly  less. 
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2.4  Linear  Hydraulic  Damping 

The  preceding  section  of  this  chapter  has  demonstrated  that  the 
drag  factor  Kn  can  be  linearized  in  a  manner  which  yields  slightly  conserva- 
tive values  of  |Z|    over  a  wide  range  of  frequencies  in  the  response 
max 

spectrum.  The  linearized  drag  factor,  IC,,  can  be  converted  into  a  linear 
hydraulic  damping  ratio,  $h»  as  shown  by  Eq.  2.7,  which  is  repeated  here  for 
convenience. 

Z  +  23  wZ  +  23hooZ  +  w2Z  =  -x   -  23.wx  (2.7) 

Equation  2.7  is  seen  to  have  the  same  form  as  Eq.  1.12  with  the  exception  of 
the  coefficients  of  the  forces  on  the  right  side  of  the  equation.  Since  the 
solution  of  Eq.  1.12  was  shown  in  Chapter  1  to  be  necessary  in  order  to  use 
the  normal  mode  method  of  solution  for  the  response  of  the  multi -degree  of 
freedom  model,  the  effects  of  3h  on  the  solution  of  Eq.  2.7  will  now  be  in- 
vestigated. 

Equation  2.7  is  solved  numerically  for  |Z|    for  all  four  earth- 
n  J  ■  'max 

quakes  with  values  of  3  ranging  from  0  to  0.10  and  values  of  3h  ranging  from 
0  to  0.10.  These  ranges  of  3  and  3,  are  considered  to  be  inclusive  of  all 
reasonable  values  of  these  variables  for  the  type  of  structure  considered  in 
this  thesis.  The  results  of  the  computations  are  shown  as  response  spectra 
in  Figs.  2.21  through  2.24.  These  response  spectra  are  seen  to  have  the  same 
general  shape  as  spectra  for  an  oscillator  without  hydraulic  damping  until 
the  high  frequency  region  is  reached.  It  should  be  noted  that  a  high  fre- 
quency tail  also  appears  in  the  spectra  for  systems  with  equivalent  linear 
hydraulic  damping.  The  magnitude  of  this  tail  can  be  shown  to  approach 
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23h  . 

( — -)(x     )  from  above.  However,  since  the  equivalent  linearization  of  the 
v  a)  /v  g  max7  ^ 

nonlinear  hydrodynamic  drag  force  has  already  been  shown  to  break  down  as 

the  drag  force  becomes  more  dominant  in  the  solution  at  high  frequencies,  an 

estimate  of  |Z|    in  the  high  frequency  tail  region  of  the  spectrum  for  a 
max 

given  system  can  best  be  made  by  the  use  of  Eq.  2.4. 

Newmark  et  al.11,ls,l%  construct  an  idealized  representation  of 
the  single  degree  of  freedom  response  spectrum  for  an  oscillator  without  drag 
forces  by  using  amplication  factors  for  displacement,  velocity,  and  acceler- 
ation. These  amplication  factors  are  multiplied  by  the  corresponding  maximum 
ground  motion  quantity  in  order  to  obtain  an  idealized  spectrum  for  an  earth- 
quake motion  of  known  (or  assumed)  values  of  x     ,  x     ,  and  x    .  If 
^  9  max   g  max      g  max 

amplification  factors  can  be  developed  for  the  response  spectrum  of  the  system 
with  equivalent  linear  hydraulic  damping  when  3  and  3h  are  known,  and  the 
high  frequency  tail  of  the  spectrum  for  the  actual  nonlinear  system  is  esti- 
mated from  Eq.  2.4,  then  an  idealized  response  spectrum  can  be  constructed  for 
a  given  set  of  conditions  for  the  oscillator  with  nonlinear  hydrodynamic  drag 
forces. 

During  the  course  of  establishing  the  amplification  factors  necessary 
for  the  construction  of  idealized  response  spectra  for  systems  with  linear 
hydraulic  damping  ratios  (Eq.  2.7),  it  was  found  that,  within  the  range  of 
variables  studied,  the  amount  of  total  linear  damping,  3  +  3.,  is  the  impor- 
tant factor  controlling  the  response  quantity  |Z|    throughout  the  spectrum 

max 

except  for  the  high  frequency  region.  Comparisons  of  |Z|    computed  from 
Eq.  2.7  at  31  frequencies  from  0.04  cps  to  15  cps  were  made  for  all  four 
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earthquake  records  for  total  damping  ratios  equal  to  .02,  .05,  .07,  .10,  .12, 

and  .15  with  the  damping  being  distributed  between  structural  and  hydraulic 

as  shown  in  Table  2.3.  The  comparison  showed  that,  up  to  the  beginning  of 

the  high  frequency  tail  in  the  spectrum,  for  a  given  value  of  the  total  damping 

ratio,  |Z|    computed  for  the  first  distribution  of  3  and  3U  was  within  5 
1  'max  s     h 

percent  of  |Z|    computed  for  the  second  distribution  of  3  and  3U-  Further, 
r        '  'max  s     n 

the  great  majority  of  the  comparisons  revealed  a  difference  of  less  than  2 
percent  between  |Z|    computed  for  the  two  different  distributions  of  damping 
at  a  given  value  of  the  total  damping  ratio.  The  ratio 

izU for  8s =  A-  Bh  ■ B 


Zl    for  3  =  B,  3,  =  A 
'max     s      h 


was  found  to  be  sometimes  greater  than  1  and  sometimes  less  than  1,  indicating 
that  at  least  a  small  part  of  the  difference  may  be  the  result  of  the  conver- 
gence criterion  (0.1  percent)  used  in  the  numerical  integration  procedure. 

A  comparison  of  |Z|  ,  when  3  =0.00  and  3U  =  0.15  with  |Z|    when 
K        '  'max      s  h  '  'max 

3  =  0.15  and  3h  =  0.00  was  made  for  each  of  the  four  earthquakes  considered. 
It  was  found  that,  for  frequencies  up  to  the  beginning  of  the  high  frequency 
tail  of  the  spectrum,  the  value  of  |Z|    computed  for  the  first  distribution 
of  total  damping  was  not  always  within  5  percent  of  the  value  of  |Z|    com- 
puted for  the  second  distribution  of  total  damping.  Most  of  the  comparisons 
differed  by  less  than  10  percent,  but  a  few  showed  differences  on  the  order 
of  15  percent. 

The  conclusion  to  be  drawn  is  that,  for  frequencies  up  to  the  be- 
ginning of  the  high  frequency  tail  region,  |Z|    computed  from  Eq.  2.7  is 

max 
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dependent  upon  the  total  linear  damping  in  the  system,  and  not  the  distribution 
of  damping  between  structural  damping  and  equivalent  linear  hydraulic  damping 
for  values  of  3h  up  to  and  including  0.10.  Thus,  the  use  of  the  amplification 
factors  published  in  Ref.  14,  and  included  in  this  thesis  as  Table  2.4,  is  a 
valid  and  proper  method  of  constructing  a  portion  of  an  idealized  response 
spectrum  for  the  single  degree  of  freedom  oscillator  with  drag  characteristics. 

The  scale  of  the  response  spectra  shown  in  Figs.  2.21  through  2.24 
is  so  small  that  the  conclusion  is  difficult  to  see  from  the  figures.  For 
illustrative  purposes,  a  comparison  of  response  spectra  for  0.  =  .07  is  shown 
in  Fig.  2.25  for  the  Ferndale  earthquake.  It  is  readily  seen  from  this  figure 
that,  up  until  the  high  frequency  tail  begins  to  appear,  the  spectral  values 
are  a  function  of  the  total  amount  of  linear  damping  in  the  system,  and  the 
relative  distribution  of  the  total  damping  between  structural  and  hydraulic 
has  little  effect. 

2.5  Idealized  Elastic  Deformation  Response  Spectrum 

It  was  mentioned  in  Section  1.5.1  that  the  characteristic  shape  of 
response  spectra  for  earthquake  motions  has  led  to  the  development  of  amplifi- 
cation factors  from  which  an  idealized  response  spectrum  can  be  constructed 
for  a  given  set  of  conditions.  The  idealized  spectrum  permits  an  estimate 
to  be  made  of  the  maximum  response  of  the  system.  This  response  quantity 
may  then  be  used  as  a  modal  response  quantity  for  a  particular  mode  of  a 
multi -degree  of  freedom  system,  and  the  various  modal  response  quantities 
may  then  be  combined  in  some  manner  to  yield  an  approximate  solution  for  the 
multi -degree  of  freedom  system  as  outlined  in  Chapter  1. 

The  principles  and  concepts  developed  in  Sections  2.2  through  2.4 
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are  summarized  here  for  the  reader  in  such  a  manner  as  to  demonstrate  the 
construction  of  an  idealized  response  spectrum  for  a  single  degree  of  freedom 
system,  with  nonlinear  hydrodynamic  drag  forces,  subjected  to  a  strong-motion 
earthquake  of  the  type  used  in  this  thesis. 

To  construct  an  idealized  deformation  response  spectrum  for  a  given 


system: 


1.  For  the  system,  compute  the  drag  factor,  IC,  and  the  undamped 
circular  natural  frequency  of  the  system,  u). 

2.  Establish  the  structural  damping  ratio,  3  ,  to  be  used  in  the 
system.  For  the  type  of  structure  discussed  in  this  thesis, 

a  structural  damping  ratio  of  between  0.005  (when  the  stress 
is  below  one-quarter  of  the  yield  stress)  and  0.05  (when  the 
stress  is  at,  or  just  below,  the  yield  stress)  is  reasonable 
to  assume.11* 

3.  Select  the  values  of  x   „  ,  x     ,  x   ,  for  which  the 

g  max   g  max   g  max 

spectrum  is  to  be  constructed. 

4.  Compute  the  quantity 

<x Jx  l> 
K  (   9  9   ) 

D  (   <x2>   ^ 

*h  "  9 

n        2  M  03 

It  should  be  noted  from  Table  1.1  that  the  magnitude  of 

^•2.:  i 

<X   X   > 

is  approximately  one-half  the  magnitude  of  x     for  the  four 

g  max 

earthquakes  shown. 
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5.  Enter  Table  2.4  with  the  value  of  6  expressed  as  a  percent, 

obtain  the  displacement  amplification  factor,  multiply  the 

displacement  amplification  factor  times  x   „  ,  and  plot  the 
r         r  g  max     r 

result  on  tripartite  logarithmic  graph  paper  as  a  displacement 
line.  This  is  shown  as  line  1  on  Fig.  2.26.  The  method  of 
linearizing  the  drag  force  by  using  ground  velocity  character- 
istics has  been  shown  to  be  unconservative  in  the  low  frequency 
region.  The  use  of  3  alone  means  that  any  hydraulic  damping 
effects  in  the  low  frequency  region  are  neglected.  This  is 
reasonable,  however,  because  the  small  amount  of  velocity  re- 
sponse in  the  low  frequency  region  results  in  minor  amounts  of 
hydraulic  damping.  Further,  as  the  frequency  of  the  system 
becomes  lower  and  lower,  the  amount  of  damping  of  any  type  has 
little  effect  on  the  maximum  deformation. 

6.  Enter  Table  2.4  with  the  amount  of  total  damping  (3  +  3.  as 
computed  above)  and  obtain  velocity  and  acceleration  amplifi- 
cation factors.  Multiply  the  velocity  amplification  factor 
times  x     and  plot  the  result  as  a  horizontal  line.  Multiply 

the  acceleration  amplification  factor  times  x     and  plot  the 

g  max    r 

result  as  an  acceleration  line.  These  two  lines  are  shown  as 
lines  2  and  3  on  Fig.  2.26. 

7.  Plot  the  quantity  (~j<-)(x  max)  as  a  displacement  line.  This 
line  is  shown  as  line  4  on  Fig.  2.26. 

8.  The  envelope  formed  by  the  intersection  of  the  four  lines  plotted 
is  the  idealized  response  spectrum  for  the  system  and  conditions 
given. 
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With  the  idealized  response  spectrum  constructed  as  described  in 
(1)  through  (8)  above,  an  estimate  of  the  maximum  response  of  the  system  is 

obtained  by  entering  the  figure  with  the  natural  frequency  of  the  system  and 

2 
reading  |Z|   ,  pseudo-velocity  (u)|Z|   ),  or  pseudo-acceleration  (w  lzlmax)« 

If  the  characteristics  of  the  system  are  changed,  a  new  response  spectrum 

should  be  constructed. 

Linear  interpolation  can  be  used  to  obtain  amplification  factors 

for  values  of  damping  not  shown  exactly  in  the  table. 

2.6  Hydrodynamic  Parameters 

The  hydrodynamic  force  exerted  by  the  fluid  on  the  mass  of  the  oscil- 
lator as  it  moves  through  the  fluid  is  represented  by  Mori  son's  formula,  Eq.  1.2 
as  described  in  Section  1.2.  This  force  expression  contains  two  dimensionless, 
empirical  coefficients,  Cx  and  Cn,  which  are  dependent  upon  the  response  of 
the  oscillator  for  their  magnitude.  Many  investigators  have  undertaken 
studies  to  determine  the  value  of  these  coefficients  under  various  conditions. 
The  results  of  some  of  these  investigations  are  published  and  discussed  in 
Refs.  17  and  18.  Morison's  investigations2  show  a  variation  of  the  drag  co- 
efficient Cn  with  Reynolds  Number,  which  is  related  to  the  relative  velocity 
between  the  mass  and  the  fluid.  Keulegan  and  Carpenter19  have  shown  a  variation 
in  both  coefficients  with  a  dimensionless  period  parameter  which  can  be  in- 
terpreted as  the  ratio  of  the  relative  displacement  between  the  mass  and  the 
fluid  to  the  pile  diameter  for  forces  against  circular  piles.  Laird20  has 
shown  that  the  two  coefficients  may  vary  greatly  for  flexible  oscillating 
cylinders.  Evaluation  of  the  two  coefficients  for  oscillatory  motion  is  further 
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complicated  by  turbulence  associated  with  eddy  formations  on  either  side  of 
the  mass  as  it  moves  back  and  forth  through  the  fluid.18  This  thesis  con- 


i 


siders  that  the  drag  factor,  Kn,  and  the  added  mass  of  the  oscillator,  m  , 
remain  constant  throughout  the  entire  time  history  of  the  oscillator  motion. 
The  values  of  IC  and  m  used  in  the  calculations  in  this  thesis  do  not  pre- 
suppose any  specific  values  for  CT  and  C~,  except  as  values  were  assigned 
to  these  coefficients  by  Burke4  in  the  formulation  of  properties  of  a  series 
of  offshore  towers  which  are  used  as  sample  structures  in  a  later  chapter. 
However,  the  designer  of  an  actual  offshore  structure,  whose  physical  pro- 
perties and  dimensions  will  eventually  be  known,  must  have  a  means  of  esti- 
mating the  coefficients  C,  and  Cn  with  reasonable  accuracy  in  order  to  analyze 
the  response  of  his  structure.  Further,  investigators  who  are  engaged  in 
experimental  research  on  the  response  of  submerged  or  partially  submerged 
bodies  to  earthquake  excitation  should  have  a  method  of  estimating  the  coef- 
ficients in  order  that  the  experimental  work  may  be  conducted  over  a  reason- 
able range  of  variables.  Since  C,  and  Cn  are  most  commonly  considered  to  be 
dependent  upon  the  relative  motion  between  the  mass  and  the  fluid,  x  and  x 
in  the  case  of  the  mass  moving  through  the  fluid  at  rest,  a  method  is  developed 

in  this  section  for  estimating  the  response  quantities  |x     ,  <|x  |>,|x     , 

r    ^        '  s'max   i  s1   '  s'max 

<|x  |>  (hereinafter  called  hydrodynamic  parameters)  for  strong  motion  earth- 
quakes of  the  type  considered  in  this  thesis. 
Consider  the  equations 

MZ  +  CZ  +  Kn.Z  +  KZ  =  -Mx   -  Km  x  (2.6) 

s       UL  g      DL  g 

and 
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Z  +  23  oi  +  2B.coZ  +  to  Z  =  -x„  -  23.wxo         (2.7) 
s       h  g      n  g 

is  which  the  nonlinear  drag  term  K-.  has  been  linearized  in  the  manner  dis- 
cussed  previously.  In  the  response  computation,  the  quantities  xg  and  xg 
may  be  computed  at  each  time  point  from  the  relations 


xc  =  Z  +  xn 
s  9 


and 


xc  =  Z  +  xn 
s         g 


Equation  2.7  is  solved  numerically  for  values  0f  3  and  3.  which  are  reasonable 

to  expect  for  the  type  of  system  considered,  and  again  it  is  found  that  the 

hydrodynamic  parameters  |x  !ma  .  <|x  l>,  |x  lav,  and  <|x  l>  are  affected  by 

s  max     s      s  max        s 

the  total  amount  of  damping  in  the  system  and  not  by  the  distribution  of  damp- 
ing between  structural  and  hydraulic.  Comparisons  were  made  of  the  hydro- 
dynamic  parameters  when  the  total  amount  of  damping  was  equal  to  2,  5,  7, 
10,  12  and  15  percent  of  critical  with  the  damping  being  distributed  between 
3  and  3.  as  shown  in  Table  2.3.  The  difference  in  a  hydrodynamic  parameter 
at  a  given  value  of  total  damping  ratio  is  within  two  to  five  percent  for  the 
different  distributions  of  damping  between  3  and  3h>  thus  indicating  that, 
if  the  amount  of  total  damping  can  be  estimated,  the  parameter  can  be  esti- 
mated. The  comparison  holds  true  for  all  four  hydrodynamic  parameters. 

Spectra  for  the  hydrodynamic  parameters  Ixl   ,  <|x  |>,  |x  I   , 
r  j         j  r  i  simax   I  si   I  simax 

and  <|x  |>  are  shown  in  Figs.  2.27  through  2.30  for  total  damping  ratios  up 
to  ten  percent.  The  general  shape  of  the  spectra  indicates  that  an  idealized 
spectrum  may  be  able  to  be  constructed  for  a  given  set  of  conditions  in  much 
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the  same  manner  as  an  idealized  deformation  response  spectrum  is  constructed. 
Each  hydrodynamic  parameter  spectrum  was  plotted,  a  smooth  spectrum  which 
approximated  the  actual  spectrum  was  established  by  the  author,  and  appro- 
priate amplification  factors  were  computed  for  each  of  the  four  earthquakes 
considered.  These  amplification  factors  were  then  averaged  in  a  mean  square 
sense  to  yield  the  amplification  factors  listed  in  Table  2.5.  The  use  of 
these  amplification  factors  to  construct  idealized  spectra  for  the  hydro- 
dynamic  parameters  |x     ,  <|x  |>,  |x  I   ,  and  <|x  |>  is  explained  in  the 
J  r  '  s'max   '  s1   '  s'max      '  s1 

next  section  of  this  chapter. 

2.7  Idealized  Elastic  Hydrodynamic  Parameter  Response  Spectra 

The  construction  of  idealized  response  spectra  for  the  hydro- 
dynamic  parameters  discussed  in  Section  2.6  is  very  similar  to  the  construc- 
tion of  an  idealized  elastic  deformation  response  spectrum. 

2.7.1  Displacement  Parameter 

To  construct  an  approximate  response  spectrum  for  |x  |    (or  <|x  |>) 
for  a  given  system: 

1.  Determine  the  quantities  Kn,  w,  3  ,   and  3h  in  the  same  manner 
as  described  in  Section  2.5. 

2.  Enter  the  appropriate  section  of  Table  2.5  with  the  value  of 
3  expressed  as  a  percent  and  obtain  the  displacement  ampli- 
fication factor.  Multiply  this  amplification  factor  times  x 

g  max 

(or  <|x  |>)  and  plot  the  result  as  a  displacement  line  on  tri- 
partite logarithmic  paper.  This  line  is  shown  as  line  1  on 
Fig.  2.31. 
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3.  Enter  the  appropriate  section  of  Table  2.5  with  the  amount  of 
total  damping  (3  +  3h)  and  obtain  the  velocity  amplification 
factor.  Multiply  this  amplification  factor  times  x     (or 
<|x  |>)  and  plot  the  result  as  a  velocity  line  as  shown  by 
line  2  on  Fig.  2.31 . 

4.  Plot  the  value  x     (or  <|x  |>)  as  a  displacement  line.  This 

g  max      g 

is  shown  as  line  3  on  Fig.  2.31. 

5.  The  envelope  formed  by  the  intersection  of  these  three  lines 

is  an  idealized  spectrum  for  |x  I    (or  <|x  |>). 
r         '  s'max     '  s1 

2.7.2  Velocity  Parameter 

To  construct  an  idealized  response  spectrum  for  |x  I  „  (or  <|x  |>) 

1  s'max     '  s1 

for  a  given  system: 

1.  Obtain  and  use  displacement  and  velocity  amplification  factors 
from  the  appropriate  section  of  Table  2.5  exactly  as  described 
above  for  the  construction  of  an  idealized  spectrum  for  a  dis- 
placement hydrodynamic  parameter.  The  resulting  lines  are  shown 
as  lines  1  and  2  on  Fig.  2.32. 

2.  Enter  the  appropriate  section  of  Table  2.5  with  the  amount  of 
total  damping  (3  +  3h)  and  obtain  the  acceleration  amplifica- 
tion factor.  Multiply  this  amplification  factor  times  x 

g  max 

(or  < | x  | > )  and  plot  the  result  as  an  acceleration  line.  This 
is  shown  as  line  3  on  Fig.  2.32. 

3.  Plot  the  value  of  x     (or  <|x  |>)  as  a  velocity  line  as  shown 
by  line  4  on  Fig.  2.32. 
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4.  The  envelope  formed  by  the  intersection  of  the  four  lines  is 

an  idealized  response  spectrum  for  |x  I  ,v  (or  <|x  |>). 
r     r         '  s'max     '  s' 

Linear  interpolation  can  be  used  to  obtain  amplification  factors 
for  values  of  damping  not  shown  exactly  in  the  table. 

Unfortunately,  no  reasonably  consistent  ratios  were  found  to  exist 

between  <|x  |>,  <|x  |>,  <|x  |>  and  other  ground  motion  quantities  such  as 

x    ,  x     ,  x     for  the  four  earthquakes  considered.  Thus,  the  de- 
g  max   g  max   g  max  ^ 

signer  of  an  offshore  structure  who  wishes  to  estimate  the  hydrodynamic  param- 
eters <|x  |>  and  <|x  |>  will  have  to  obtain  the  values  of  <|x  |>,  <|x  |>,  and 

I   Sl  I   Sl  I   gl        I   gl 

<|x  |>  from  processing  an  actual  earthquake  record  or  estimating  them  using 

the  values  shown  in  Table  1.1  as  a  guide.  The  range  of  the  ratios  between 

<|x  |>,  <|x  |>,  <|x  |>  and  x     ,  x     ,  x     are  repeated  here  for  the 
1  g1     g      g       g  max   g  max   g  max      v 

convenience  of  the  reader. 

Ratio  Lowest  Value        Highest  Value 


0.194  0.385 


0.124  0.167 


0.084  0.185 


It  should  be  noted  that  the  idealized  spectra  for  hydrodynamic  param- 
eters become  less  descriptive  of  the  oscillator  response  at  very   low  frequencies, 
about  0.08  cps  and  less  for  the  four  earthquakes  used  in  the  computations.  This 

is  to  be  expected,  since  the  quantities  |x     »  <|x  |>,  |x  I   ,  and  <|x  |> 

1  s'max   '  s1   '  s'max      '  s1 

become  quite  small  for  very   soft  (low  frequency)  systems. 
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2.8  Inelastic  Response  Spectra 

In  this  section,  the  effects  of  the  nonlinear  hydrodynamic  drag 
force  on  the  response  of  the  oscillator  in  the  inelastic  range  are  investi- 
gated. Figure  2.33  shows  the  resistance-deformation  relationship  of  the 
spring  in  a  system  with  elasto-plastic  spring  characteristics.  The  character- 
istics of  the  elasto-plastic  system  are  an  initial  elastic  region,  a  constant 
plastic  resistance  limit,  and  an  elastic  unloading  region.  The  unloading 
portion  is  considered  to  have  the  same  slope  as  the  initial  elastic  curve,  and 
the  unloading  is  elastic  until  yielding  is  reached  in  the  opposite  direction.12 
The  yield  point  deformation  is  denoted  by  Z  ,  and  the  absolute  maximum  dis- 
placement of  the  mass  relative  to  the  base  is  denoted  by  |Z|  ,  .  The  yield 

1  'max 

resistance  is  denoted  by  R  .  The  ductility  factor,  y,  is  defined  as 

y  =    zmax  (2.9) 

y 

and  expresses  the  maximum  deformation  of  the  inelastic  system  in  terms  of  its 
yield  point  deformation.  With  y  defined  in  this  manner,  the  maximum  deforma- 
tion of  the  system  is  yZ  ,  the  maximum  elastic  deformation  is  Z  ,  and  the 
maximum  inelastic  deformation  is  (y  -  1)Z  .  A  variety  of  calculations  of 
the  response  of  inelastic  systems  has  been  made  and  reported  in  the  litera- 
ture. 11_ltt  The  results  of  the  calculations  are  conveniently  displayed  in 
the  form  of  response  spectra  in  much  the  same  manner  as  response  spectra 
for  elastic  systems  are  displayed.  For  elasto-plastic  systems,  response 
spectra  can  be  plotted  for  constant  values  of  the  ductility  factor,  y,  and 
when  the  plot  is  made  in  terms  of  the  elastic  component  of  the  displacement, 
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the  accelerations  are  properly  presented.12  In  this  type  of  spectrum  represen- 
tation, the  pseudo-velocity  quantity  plotted  is  wZ  .  Elasto-plastic  response 
spectra  for  systems  with  3C  =  0.02,  Knx  __Y/K  =  0.00  (no  drag),  subjected 
to  the  El  Centro  earthquake  are  shown  in  Fig.  2.34  for  values  of  u  =  1  (elastic 
case)  to  u  =  4.  In  these  spectra,  the  total  deformation  of  the  system  can  be 
determined  by  reading  the  displacement  quantity,  which  is  Z  ,  and  then  using 
the  relationship 

'  'max    y  y 

In  the  low  frequency  range  of  the  spectra  shown  in  Fig.  2.34,  the 
total  deformation  is  approximately  the  same  for  the  elasto-plastic  systems 
as  for  the  elastic  systems.  In  the  high  frequency  region,  the  pseudo-acceler- 
ation is  approximately  the  same  for  the  elasto-plastic  systems  as  for  the 
elastic  system.  In  the  mid-frequency  region,  some  generalizations  have  been 
made  regarding  the  amount  of  reduction  in  yield  level  for  different  values 
of  y,  but  more  detailed  generalizations  are  still  under  study.12 

Elasto-plastic  deformation  response  spectra  for  systems  with  non- 
linear hydrodynamic  drag  forces  subjected  to  the  El  Centro  earthquake  are 
shown  in  Figs.  2.35  through  2.37.  Again,  in  the  low  frequency  region,  the 
total  deformation  is  approximately  the  same  for  the  elasto-plastic  systems 
as  for  the  elastic  system.  In  the  mid-frequency  region,  the  spectra  have 
the  same  general  shape  as  do  the  elasto-plastic  spectra  for  systems  without 
the  drag  forces,  but  no  quantification  of  the  results  has  been  made.  Further 
study  in  the  direction  of  establishing  rules  for  the  construction  of  idealized 
elasto-plastic  response  spectra  in  this  frequency  region  is  clearly  indicated. 
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In  the  high  frequency  region,  the  elasto-plastic  spectra  for  the  systems  with 
the  lowest  value  of  the  drag  parameter  (0.006)  approach  approximately  the  same 
values  in  the  high  frequency  tail  as  the  elastic  spectrum  does.  However,  for 
the  systems  with  the  greatest  value  of  the  drag  parameter  (0.50),  the  high 
frequency  tails  of  the  elasto-plastic  spectra  approach  values  which  are  re- 
duced by  as  much  as  50  percent  from  the  elastic  case.  In  other  words,  there 
is  a  reduction  in  the  magnitude  of  the  elasto-plastic  high  frequency  tail  from 
that  obtained  in  the  elastic  case,  and  the  amount  of  the  reduction  is  a  func- 
tion of  the  amount  of  nonlinear  hydrodynamic  drag  in  the  system.  In  the  sys- 
tems with  a  small  amount  of  drag,  the  reduction  is  yery   slight,  and  may  not 
merit  much  consideration,  but  the  reduction  for  high  drag  systems  can  be  sig- 
nificant. A  study  of  the  numerical  results  obtained  in  this  investigation 

revealed  that,  for  the  spectra  with  3  =  0.02,  Knx    /K  =  0.50,  values  of 

K  s        D  g  max 

u  up  to  4  and  the  El  Centro  earthquake,  the  factor  ~ — ^-  applied  to  the 
elastic  case  high  frequency  tail  yields  a  satisfactory  estimate  of  the  elasto- 
plastic  high  frequency  tail  (see  Table  2.6).  The  reader  is  cautioned  that 
this  expression  for  the  reduction  factor  is  applicable  only  to  the  specific 
systems  and  earthquake  used,  and  is  not  to  be  extended  to  other  conditions. 

In  order  to  compare  the  effects  of  the  nonlinear  hydrodynamic  drag 
force  on  the  elasto-plastic  response  spectra,  spectra  are  plotted  for  given 
values  of  structural  damping  ratio,  3  ,  and  ductility  factor,  u,  for  systems 
with  varying  amounts  of  drag  subjected  to  the  El  Centro  earthquake.  These 
comparisons  are  shown  in  Figs.  2.38  through  2.40.  These  figures  show  that, 
for  a  small  ductility  factor  (u  =  1.5),  increasing  amounts  of  nonlinear  hydro- 
dynamic  drag  reduce  the  value  of  the  yield  deformation,  Z  ,  for  a  range  of 
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frequencies  between  low  frequencies  and  the  beginning  of  the  high  frequency 
tail.  This  means  that  the  maximum  response  of  the  system  with  high  drag  forces 
is  less  than  the  maximum  response  of  the  system  with  lower  drag  forces,  or, 
the  drag  forces  tend  to  produce  a  damping  effect  on  the  maximum  deformation 
of  the  spring  in  the  mid-frequency  region  of  the  spectrum.  For  example,  in 
Fig.  2.38,  at  a  frequency  of  1  cps,  and  the  highest  value  of  the  dimensionless 
drag  parameter  used  (0.50),  Z  is  only  about  50  percent  of  the  value  of  Z 
for  the  system  with  no  drag.  However,  as  the  ductility  factor  u  increases, 
the  effects  of  the  nonlinear  hydrodynamic  drag  force  become  less  and  less 
noticeable  in  the  mid-frequency  range.  In  Fig.  2.40,  with  u  =  4,  the  value 
of  Z  for  the  highest  drag  system  is  about  90  percent  of  the  value  of  Z  for 
the  system  with  no  drag  at  a  frequency  of  1  cps.  This  is  reasonable  to  expect, 
since  the  velocity  of  the  oscillator,  x  ,  can  be  expected  to  be  smaller  during 
plastic  deformation  of  the  spring  than  during  elastic  deformation  of  the 
spring,  and  the  magnitude  of  the  drag  force  is  a  function  of  the  oscillator 
velocity.  It  is  also  to  be  noted  from  Fig.  2.40  that  it  is  possible  for  Z 
in  the  case  of  the  system  with  drag  to  exceed  slightly  the  value  of  Z  for 
the  system  with  no  drag  at  the  same  frequency  (note  the  response  values  at 
frequencies  of  0.8  cps  and  3.0  cps)  in  the  mid-frequency  range.  This  can  re- 
sult since  the  maximum  response  of  the  elasto-plastic  system  depends  on  the 
entire  history  of  the  deformation  and  the  energy  absorbed  in  the  spring  during 
plastic  deformation  as  well  as  the  energy  dissipated  by  the  damping  effects 
of  the  fluid  surrounding  the  oscillator. 

The  significant  features  of  the  elasto-plastic  response  spectra 
for  systems  with  nonlinear  hydrodynamic  drag  are: 
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1.  The  drag  does  not  have  much  effect  on  the  spectrum  in  the  low 
frequency  region; 

2.  A  high  frequency  tail  does  exist,  and  large  amounts  of  drag 
can  cause  significant  reductions  in  Z  for  increasing  values 

y 

of  the  ductility  factor  (less  significant  reductions  in  Z  with 
increasing  values  of  u  for  low  drag  systems); 

3.  The  damping  effects  of  the  drag  in  the  mid-frequency  region 

of  the  spectrum  become  less  pronounced  as  the  ductility  factor 
increases,  even  for  the  systems  with  high  drag. 
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CHAPTER  3 
MULTI-DEGREE  OF  FREEDOM  SYSTEMS 

3.1  General 

In  this  chapter,  the  principles  developed  and  discussed  in  Chapter 
2  for  the  single  degree  of  freedom  oscillator  are  extended  to  the  multi- 
degree  of  freedom  system.  The  nonlinear  drag  forces  are  linearized  in  such 
a  manner  as  to  insure  the  existence  of  classical  normal  modes  in  the  equa- 
tions of  motion  in  order  that  normal  mode  superposition  may  be  used  as  an 
approximate  solution  to  the  problem.  The  method  of  linearization  is  shown 
to  be  accurate  for  a  series  of  five  degree  of  freedom  systems  which  are 
representative  of  the  type  of  structure  being  investigated.  A  means  of  es- 
timating the  amount  of  equivalent  linear  hydraulic  damping,  using  ground 
motion  velocity  characteristics,  is  developed  and  is  shown  to  yield  good  re- 
sults. The  construction  of  an  idealized  response  spectrum  for  each  mode  of 
the  multi -degree  of  freedom  system  is  discussed  and  described,  and  sample 
problems  are  solved  by  the  superposition  of  normal  mode  responses. 

3.2  Approximate  Solution  of  the  Nonlinear  Equations  of  Motion 

In  this  section,  an  approximate  method  of  solving  the  nonlinear 
equations  of  motion  for  the  multi-degree  of  freedom  system  is  developed  and 
shown  to  yield  satisfactory  results.  Various  features  of  the  solution,  such 
as  the  formulation  of  the  structural  damping  matrix,  the  linearization  of 
the  drag  forces,  the  superposition  of  modal  responses,  and  the  response  of 
very   high  frequency  modes,  are  discussed.  The  transformation  of  the  nonlinear 
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equations  of  motion  into  uncoupled  linear  modal  equations,  which  was  first 
shown  in  Chapter  1,  is  repeated  in  this  section  for  the  convenience  of  the 
reader. 

3.2.1  Equations  of  Motion 

The  equations  of  motion  for  the  lumped  parameter  model  shown  in 
Fig.  1.1(b)  are  written  in  matrix  form  as  Eq.  1.3  and  are  repeated  here  to 
form  the  basis  of  the  analytical  discussion  in  this  section. 

FMJ{Z}  +  [Cs]{Z}  +  [K]{Z}  =  -r;Mj{xg}  -  rKD^  {xs|xs|}   (3.1) 

No  structural  dampers  are  shown  on  the  model,  but  the  damping  matrix  [C  ] 
can  take  a  variety  of  forms  as  will  be  discussed  in  later  paragraphs.  For 
a  given  structural  system  and  known  earthquake  record,  response  quantities 
of  interest  may  be  obtained  by  the  numerical  integration  of  Eq.  3.1,  but 
this  procedure  is  costly  and  requires  the  use  of  a  digital  computer.  If 
Eq.  3.1  can  be  satisfactorily  approximated  by  a  system  of  simultaneous 
linear  differential  equations  of  motion  which  possess  classical  normal 
modes,  then  a  solution  can  be  obtained  by  the  superposition  of  normal  mode 
solutions.  A  further  approximation  can  be  made  by  using  response  spectra 
to  obtain  the  maximum  response  in  each  mode  and  then  combining  the  modal 
maxima  in  some  fashion,  s/ery   much  like  the  approximate  modal  techniques, 
using  response  spectra,  which  are  reported  in  the  literature11"13  for  land- 
based  structures. 

In  the  normal  mode  method,  the  vector  of  a  response  quantity  ex- 
pressed in  one  coordinate  system,  such  as  the  relative  displacement  between 
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each  node  and  the  base,  {Z},  is  transformed  into  the  "normal"  coordinate  sys- 
tem by  the  relation 

{Z}  =  [<j>]{Y}  (3.2) 

If  the  nonlinear  drag  force  were  not  present,  Eq.  3.2  and  the  ortho- 
gonality relationship  between  [<J>]  and  ["Mj,  [K],  and  [C  ]  could  be  used  to 
transform  Eq.  3.1  into  nm  (number  of  modes)  uncoupled  single  degree  of  free- 
dom equations  in  the  normal  coordinate  system.  If  the  nonlinear  drag  force 
can  be  converted  into  an  equivalent  linear  drag  force,  then  the  resulting 
linear  differential  equations  of  motion  can  be  transformed  into  uncoupled 
single  degree  of  freedom  equations  of  motion  in  the  normal  coordinate  system 
as  follows:  Consider  that  ["Kn-J  ^xslx  U  can  De  transformed  into  an  equiva- 
lent linear  force,  IT" Kni  ^  {**}     i.e., 

n<D^]{xs|xs|}  -   rKD|_J  (xs>  (3.3) 

Substitution  of  Eq.  3.3  into  Eq.  3.1  yields, 

FMJ{Z}  +[CS]{Z}+  [K]{Z}  =  -FMj{xg}  -  [-KDLJ{xs}    (3.4) 
Since  {Z}  =  {x  }  -  {x  },  Eq.  3.4  can  be  rewritten  as 

[^MJ{Z}  +  [CS]{Z}  +  rKDL-J{Z}  +  [K]{Z} 

=  -FMj{xg}  -  rKD1/J{xg}  (3.5) 

Substitute  Eq.  3.2  into  Eq.  3.5  and  premultiply  the  result  by  M 
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MVmJEcDKY}  +  [cj)]T[Cs][(()]{Y}  +  [<f>]TrKDLJ|>]{Y}  +  [<f>]T[K]|>]{Y} 
=     -U]VMj{xg}  -   [<f)]TrKDL-J{xg}  (3.6) 

Recall   the  following  definitions  from  Chapter  1: 

*  T 

C"M  -J     =     [<J>]   M[(J>]  =     a  diagonal  modal  mass  matrix 

*  T 

[~~C  -J  =  M   [C  ~\M         =  a  diagonal  modal  structural  damping 

matrix.  The  form  of  [C  ]  which  permits 

this  is  discussed  later. 

*  T 

["K  ~J     =     [<|>]   [K][<|>]  =     a  diagonal  modal   stiffness  matrix 

%  T 

[Kn(  ]       =     [<j>]   C~~ Kpj.  -J[<{>]  =     a  modal   linear  hydraulic  damping 

matrix  which  is  not  necessarily 
diagonal 
C"oo  J     =     [*M  ^]     ["K  ~J  =     a  matrix  of  circular  natural   fre- 
quencies of  the  system 

Substitution  of  these  definitions  into  Eq.   3.6  yields 

FM*J{Y}  +  rC*^]{Y}  +  [KpL]{Y}  +  rK*J{Y} 

=     -Dofr-MjCXg}  -  [(0]TrKDLJ{xg}  (3.7) 

At  this  point,  the  formulation  of  the  structural  damping  matrix,  [C  ],  and 
the  linearization  of  the  drag  factor,  r~Kn-]>  will  be  discussed. 

3.2.2  Structural  Damping 

If  the  structural  damping  matrix  [C  ]  is  proportional  to  the  mass 
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matrix,  then  [~C  -J  will   be  diagonal   since  the  natural  mode  shapes  are  ortho- 
gonal with  the  mass  matrix.     This  formulation  of  [Cs]  results  in  modal  damping 
ratios  which  decrease  with  increasing  natural   frequency.     This  is  shown  as 
follows: 
Let, 

[C  ]     =     CFMJ 


then, 


Define 


then, 


where 


[4)]T[Cs][d>]     =     S[<J>]T^Mj[cf>]     =     cr-M*^]     =     rc*J  (3.8) 


rM*J_Vc*^l     =     r~23so)J     =     ?[I]  (3.9) 


3        =     o*—  (3.10) 

2wii 


$    =  structural  damping  ratio  in  the  ith  mode 
si 

a>. .  =  circular  natural  frequency  of  the  ith  mode 

In  a  similar  manner,  it  can  be  shown  that,  if  the  structural  damping 

•k 

matrix  [C  ]  is  proportional  to  the  stiffness  matrix  [K],  then  [~  C  J  is  a 
diagonal  matrix  in  which  the  modal  damping  ratios  3c  increase  with  increasing 
frequency. 

The  elements  of  the  structural  damping  matrix  [C  ],  can  be  selected 
to  yield  the  same  structural  damping  ratio  in  each  mode.  It  is  common  prac- 
tice to  assume  the  same  structural  damping  ratio  in  all  modes  in  the  dynamic 
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analysis  of  damped  multi-degree  of  freedom  systems.12  If  the  modal  matrix 
[<{>]  is  normalized  such  that 

[4>]TrMJ[>]  =  [MJ     (a  unit  matrix)  (3.11) 

Then 

[Cs]  =  C-MJU]rc*J[(D]T|^Mj  (3.12) 

And,  since 

r^  =  ^M*^r2esoJJ  (3.13) 

Equation  3.12  can  be  written  as 

[Cs]  =  rMj[4>]rM*^jr23sooJ[cJ>]TrMj  (3.14) 

Thus,  it  is  seen  that  [C  ]  can  be  formulated  in  such  a  manner  as 
to  yield  modal  structural  damping  ratios  which  are  constant,  which  increase 
with  increasing  mode  number,  or  which  decrease  with  increasing  mode  number. 
In  the  computations  reported  in  this  thesis,  the  structural  damping  matrix 
[C  ]  is  computed  in  such  a  way  that  3  is  the  same  in  all  modes. 

3.2.3  Linearization  of  the  Drag  Force 

Foster,9  and  Malhotra  and  Penzien6  linearize  the  nonlinear  drag 
forces,  ["HL-J{x  |x  |}  by  replacing  them  with  a  velocity  dependent  linear 
force,  r~KD. -J{x  }  at  each  node  which  performs  the  same  amount  of  work  during 
a  period  of  time  as  does  the  nonlinear  force  at  that  node.  Foster  then 
recognizes  that  the  resulting  [""K^.-J  matrix  is  not  orthogonal  with  the  mode 
shapes  and  proceeds  to  solve  the  equations  of  motion  by  a  nonclassical  mode 
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approach  using  complex  eigenvalues  and  eigenvectors.  Malhotra  and  Penzien 
transform  ["Km-J  into  [Kn|  ]  and  then  force  [KR|  ]  to  be  diagonal  by  means 
of  an  iterative  process  that  replaces  the  nondiagonal  matrix  with  a  diagonal 
matrix  which  dissipates  the  same  amount  of  energy  over  a  period  of  time. 

The  hydraulic  drag  forces,  whether  linear  or  nonlinear,  may  be  con- 
sidered to  act  as  dampers,  that  is,  they  cause  energy  to  be  dissipated  as 
the  system  oscillates  through  the  fluid.  In  the  approximate  solution  being 
developed  in  this  section,  the  nonlinear  dampers  in  the  multi -degree  of 
freedom  system  are  replaced  with  linear  dampers  (except  for  \/ery   high  fre- 
quency modes  which  are  discussed  in  a  later  section)  which  dissipate  the 
same  amount  of  energy  over  the  total  system,  but  do  not  necessarily  preserve 
the  dissipation  of  energy  at  each  node.  Further,  the  resulting  "equivalent" 
linear  damper  at  each  node  is  made  to  be  proportional  to  the  lumped  mass  of 
that  node  in  order  that  the  modal  equations  of  motion  will  be  uncoupled.  For 
an  offshore  structure  with  a  portion  of  the  structure  above  the  water  line 
as  shown  in  Fig.  1.1(a),  this  procedure  results  in  a  linear  hydraulic  damper 
being  provided  at  the  level  of  the  top  mass  where,  in  fact,  no  nonlinear 
hydraulic  damper  really  exists.  With  this  exception,  the  proportioning  of 
the  linear  hydraulic  dampers  in  the  same  ratio  as  the  mass  is  distributed  is 
a  reasonable  approach  for  the  type  of  structure  being  discussed,  since  the 
nonlinear  damper  is  proportional  to  the  projected  area  of  the  structure  lumped 
at  the  node,  and  the  magnitude  of  the  mass  at  a  given  node  will  also  be  gen- 
erally proportional  to  the  projected  area.  The  method  of  linearizing  is  as 
follows: 
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nn 


t 


Energy  dissipated  by     =  (Kn      x     |x     |)(x     )  dt 

nonlinear  dampers  >-, 


o 


D..  s.1  s . '   s- 
JJ   J   J     3 


nn  9 

=  t   7  Kn   <)C  |x   |>                (3.15) 

•  S  D. .   s  . '  s  • ' 

J=l  JJ    J   J 


Energy  dissipated  by     nn 
equivalent  linear    =   £ 


(Knl   x  )(x  )  dt 
v  DL  . .  s .   s. 


dampers  j  =  l  J0  jj   j    j 


=  t  f     KDL   <x*>  (3.16) 

J=l     JJ    J 

where 

nn  =  number  of  nodes  in  the  system 
j   =  jth  node  in  the  system 

Equating  the  energy  dissipated, 

nn       9  nn        ? 

I       K    <kl    l*s  \y     =       I       KDL  KK   >  (3J7) 

j=l   °jj   Sj  Sj      J=l     Jj   Sj 

It  is  noted  that  only  the  total  energy  dissipated  is  preserved  and  the  dis- 
tribution of  the  energy  dissipation  among  nodes  is  not  preserved.  Now,  by 
forcing  the  linear  dampers  to  be  proportional  to  the  mass  matrix, 

rKDL^J  =  aC^Mj  (3.18) 


and  substituting  Eq.  3.18  into  Eq.  3.17 
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a  = 


nn       ~ 

J=l  JJ    J   J 
nn      .o 

y  m..  <xc  > 


(3.19) 


Recall  that  [K^]  was  defined  as  U]Vkdl~J[>].  Substituting  Eq.  3.18  into 
this  definition, 

[K^_]  =  a[<D]TrMj[>]  =  afM*.J  (3.20) 

Substituting  Eq.  3.20  into  Eq.  3.7  yields 

r-M*J{Y}  +  |^C*J{Y}  +  a^M*J{Y}  +  [~  K*J{Y} 

=  [*]T|^MJ{-xg  -  ctxg}  (3.21) 

Equation  3.21  represents  nm  single  degree  of  freedom  uncoupled  equations  of 
motion  in  the  normal  coordinate  system. 

Before  proceeding,  the  method  of  linearizing  described  by  Eqs.  3.18 
and  3.19  will  be  tested  for  accuracy  in  computing  the  response  of  multi -degree 
of  freedom  systems  to  earthquake  excitation.  A  series  of  five  degree  of  free- 
dom systems,  representative  of  offshore  towers  in  mass,  stiffness,  and  drag 
distribution,  is  developed  in  order  to  test  the  linearization  technique. 

Seven  structural  systems,  consisting  of  different  combinations  of 
mass  and  stiffness  distribution,  are  considered.  Each  system  is  further 
subdivided  into  one  or  more  distributions  of  the  drag  factor,  Kn,  at  each  node. 
The  characteristics  of  these  systems  are  shown  in  Table  3.1.  In  the  table, 
the  values  shown  for  the  mass  at  each  node  are  considered  to  be  the  total 
mass  at  that  node,  and  include  the  added  mass  term.  The  different 


60 


distributions  of  drag  factor  are  labeled  L,  M,  and  H  in  the  table,  and  the 
column  under  each  of  these  headings  represents  the  relative  magnitude  of  IC 
at  each  node  for  the  particular  distribution  shown.  For  example,  drag  factor 
distribution  M  represents  the  distribution  of  KQ  at  each  node  as  0,  1,  1,  2, 
3  at  nodes  1  through  5  respectively.  The  values  of  ZKL/Zm  shown  in  the 
table  for  the  distributions  of  drag  factor  are  of  sufficient  range  to  span 
reasonable  values  of  this  quantity  for  actual  offshore  towers  of  the  type 


being  investigated.  The  quantity  IIC/Em  for  the  four  towers  reported  by 
Burke  and  Tighe"  range  from  a  low  of  0.0034  inches"  to  a  high  of  0.0054 
inches"  .  The  relative  magnitudes  of  the  undamped  natural  frequencies  of 
the  system,  normalized  to  the  fundamental  frequency,  are   as  shown.  In  the 
computations  which  are  discussed  in  the  following  paragraphs,  one  of  the 
natural  frequencies  is  set  to  a  desired  value  and  the  fundamental  frequency 
is  computed  from  the  ratio  of  natural  frequencies.  The  value  of  k,  is  com- 
puted from  f,  and  the  quantity  f,//k,/m, ,  then  the  remaining  spring  stiffnesses 
are  computed  from  k,  and  the  distribution  of  k?  through  k^.  The  stiffness 
matrix  [K]  for  the  entire  system  is  then  assembled  and  is  a  tri -diagonal 
matrix  for  the  systems  shown. 

Each  of  the  systems  shown  in  Table  3.1  is  subjected  to  the  El  Centro 
earthquake  and  response  quantities  are   computed  by  the  numerical  integration 
of  the  equations  of  motion  using  the  criteria  for  convergence  and  time  inter- 
val described  in  Chapter  1.  For  each  system,  computations  of  maximum  re- 
sponse are  made  for  two  values  of  3  in  each  mode,  0.01  and  0.05,  and  three 
frequency  ranges.  The  low  frequency  range,  range  number  1,  is  established 

by  setting  the  fifth  natural  frequency  equal  to  the  frequency  at  which  x 

g  max 
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and  x     intersect  when  plotted  on  tripartite  logarithmic  paper.  The 
a    g  max 

mid-frequency  range,  range  number  2,  is  established  by  setting  the  funda- 
mental  frequency  equal  to  the  frequency  at  which  x     and  x     inter- 
sect. The  high  frequency  range,  range  number  3,  is  established  by  setting 
the  fundamental  frequency  equal  to  the  frequency  at  which  x     and  x 
intersect  when  plotted  on  tripartite  logarithmic  paper.  The  maximum  dis- 
placement of  each  node  relative  to  the  base,  {|Z|   },  is  computed.  The 

max 

results  of  the  computations  are  shown  in  Tables  3.2  through  3.27.  The 
column  in  the  tables  titled  "No  Drag"  shows  the  results  of  the  computations 
when  ["K^-J  is  set  to  zero.  The  column  in  the  tables  titled  "Actual  Drag" 
shows  the  results  of  the  computations  when  ["^.J  is  computed  from  the  values 
shown  in  Table  3.1,  thus  nonlinear  hydrodynamic  drag  forces,  [""Kn^]{x  |x  |}, 
are  acting  on  the  structure.  In  the  computations  for  the  "Actual  Drag" 
systems,  the  time  averages  {<x  |x  |>}  and  {<x  >}  are  collected.  These  time 
averages  are  then  used  to  linearize  the  drag  factor  and  set  it  to  be  pro- 
portional to  the  mass  matrix  as  prescribed  by  Eqs.  3.18  and  3.19.  Equation 
3.4,  using  the  linear  drag  factor  matrix  just  computed,  is  then  integrated 

numerically  to  obtain  the  response  quantity  {|Z|   }.  The  results  of  this 
*  v  M     J     '  'max 

integration  are  labeled  "Case  1"  in  the  tables  of  results.  The  column 
labeled  "Case  1/Actual  Drag"  in  Tables  3.2  through  3.27  represents  the  com- 
parison of  the  solution  of  Eq.  3.1  with  the  solution  of  Eq.  3.4.  A  study 
of  these  results  reveals  that  the  linearization  technique  described  by  Eqs. 
3.18  and  3.19  is  extremely  accurate.  The  difference  between  the  linear  and 
the  nonlinear  results  is  less  than  6  percent  (with  one  exception — system  1M, 
3S  =  0.01,  node  3,  mid-frequency  range)  and  is  less  than  2  percent  in  the 
majority  of  the  cases. 
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The  maximum  distortion  between  two  adjacent  nodes  is  also  computed 

for  both  the  linear  and  the  nonlinear  equations  of  motion.  The  results  of 

the  computations  for  the  linear  case  in  which  the  linearization  is  carried 

out  according  to  Eqs  3.18  and  3.19  are  not  shown  in  this  thesis,  but  the 

linear  solution  is  as  close  to  the  nonlinear  solution  as  for  the  response 

quantity  {|Z|   }  shown  in  the  tables.  Thus,  it  is  concluded  that  the  meth- 
max 

od  of  linearizing  the  drag  force  described  by  Eqs.  3.18  and  3.19  yields 
satisfactory  results. 

The  primary  difficulty  in  using  Eqs.  3.18  and  3.19  to  linearize  the 

•2,  • 
drag  force  lies  in  the  ability  to  determine  the  quantities  {<x  |x  |>}  and 

•  2 
{<x  >}  a  priori.  If  these  quantities  can  be  estimated  in  some  fashion,  and 

if  the  estimate  can  be  shown  to  be  accurate  in  certain  frequency  ranges, 
then  Eq.  3.1  can  be  linearized  into  Eq.  3.4,  and  Eq.  3.4  can  be  transformed 
into  uncoupled,  modal,  single-degree  of  freedom  equations  in  the  normal  co- 
ordinate system  which  can  be  solved  by  spectrum  techniques.  It  was  demon- 
strated in  Chapter  2  that,  for  the  single  degree  of  freedom  oscillator,  the 

•2. •  ,      *2 
quantities  <x  x  >  and  <x  >  are  conservative  estimates  of  the  oscillator 

g1  g1     g 

velocity  in  all  but  the  low  frequency  region  of  the  spectrum.  Since  the 
major  effects  of  hydraulic  damping  occur  in  the  velocity  amplification  re- 
gion of  the  spectrum,  it  seems  reasonable  to  estimate  the  quantities 

t  |x  |>}  and  {<x  >}  as  {<x 
quantities,  Eq.  3.19  becomes 


{<x  |x  |>}  and  {<x  >}  as  {<x  |x  |>}  and  {<x  >}  respectively.  Using  these 


a  = 


/  T     V  \/<x^|xg|>       

(3.22) 


j=l    j> 


nn 


1 
\    j=l 


•2 

i .   /  \  ^x 

jj  /  \  g 


M..  /  \  <x_  > 
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A  series  of  computations,  labeled  "Case  2"  in  Tables  3.2  through  3.27,  is 
made  in  which  the  nonlinear  drag  factor  is  linearized  in  the  manner  pre- 
scribed by  Eqs.  3.18  and  3.22.  A  study  of  the  column  labeled  "Case  2/Actual 
Drag"  shows  that  the  linearization  described  by  Eqs.  3.18  and  3.22  is,  in 
general,  conservative  in  the  mid  and  high  frequency  ranges  and  is  generally 
unconservative  in  the  low  frequency  range.  The  degree  of  conservatism  in 
the  mid-frequency  range  is  seen  to  be  less  than  5  percent  in  most  cases, 
although  one  case  is  conservative  by  12  percent.   In  the  high  frequency 
range,  the  results  are  seen  to  be  conservative  by  as  much  as  20  percent  in 
some  cases,  especially  in  the  high  drag  systems,  but  many  of  the  results  are 
conservative  by  less  than  10  percent  in  the  high  frequency  range.  The  degree 
of  unconservatism  in  the  low  frequency  range  is  seen  to  be  less  than  10  per- 
cent in  most  cases,  but  can  range  as  high  as  23  percent.  This  unconservatism 
in  the  low  frequency  range  is  to  be  expected,  since  very  soft  (low  frequency) 
systems  can  be  expected  to  exhibit  little  velocity  response  at  each  node, 
and  so  the  use  of  the  ground  velocity  characteristics  overestimates  the 
velocity  of  each  node.  One  further  observation  to  be  made  from  a  study  of 
the  numerical  results  is  that  the  use  of  the  ground  velocity  to  linearize 
the  drag  force  is  generally  more  conservative  as  the  response  is  reduced 
to  a  greater  degree  by  the  actual  nonlinear  drag  force.  In  those  cases  in 
which  the  Case  2  linearization  is  unconservative  in  frequency  ranges  2  and 
3,  the  degree  of  unconservatism  is  generally  very  slight  (less  than  1  percent). 

The  maximum  deformation  between  adjacent  nodes  is  also  computed 
and  a  comparison  of  the  linear  drag  response  (linearized  by  the  ground  velocity 
characteristics)  with  the  nonlinear  drag  response  is  shown  in  Figs.  3.1 
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through  3.13.  These  figures  show  that  the  maximum  spring  deformation  com- 
puted when  the  drag  force  is  linearized  by  the  ground  velocity  characteristics 
is  also  conservative  in  the  mid  and  high  frequency  ranges  and  generally  un- 
conservative  in  the  low  frequency  range.   It  is  concluded  from  these  compu- 
tations that  the  use  of  Eqs  3.18  and  3.22  is  a  reasonable  way  to  linearize 
the  drag  force  and  yields  slightly  conservative  results  except  in  the  low 
frequency  range.  The  conservatism  introduced  by  the  linearization  is  seen 
from  the  figures  to  be  generally  less  than  10  percent  in  the  mid  and  high 
frequency  ranges,  with  some  of  the  results  in  the  high  frequency  range  being 
conservative  by  as  much  as  20  percent.  The  results  of  the  low  frequency  range 
calculations  are  generally  unconservative  by  less  than  20  percent,  but  can  be 
unconservative  by  as  much  as  25  percent. 

The  effect  of  structural  damping  on  the  response  is  quite  apparent 
from  the  tables  and  figures.  For  example,  the  maximum  relative  displacement 
between  the  bottom  node  (mass  number  5)  and  the  ground  for  system  4M  in  fre- 
quency range  2  is  reduced  approximately  17  percent  by  the  hydraulic  drag  when 
the  structural  damping  ratio  is  0.01  in  each  mode  (see  Table  3.14).  Yet  the 
response  is  reduced  by  less  than  2  percent  by  the  hydraulic  drag  when  the 
structural  damping  ratio  is  0.05  in  each  mode  (see  Table  3.15).  Thus,  the 
selection  of  structural  damping  ratios  to  use  in  the  response  computations 
is  seen  to  be  yery   important. 

3.2.4  Modal  Equations  in  Normal  Coordinates 

The  use  of  Eqs.  3.18  and  3.22  to  linearize  the  drag  force  has  been 
shown  to  yield  generally  conservative  results  in  the  mid  and  high  frequency 
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ranges,  and  unconservative  results  in  the  low  frequency  range.  Since  the 
linearization  prescribed  by  Eqs.  3.18  and  3.22  insures  that  the  equations  of 
motion  for  the  linearized  system  (Eq.  3.5)  will  possess  classical  normal  modes, 
the  discussion  of  the  superposition  of  normal  modes  is  resumed.  Recall  the 
uncoupled  equations  of  motion  in  the  normal  coordinate  system. 

=     WTr-M^{-xg  -  axg}  (3.21) 

_     *       i 
Premulti plying  by  [/M  -~J      yields 

{Y}  +  r~23su)J{Y}  +  r26hwJ{Y)  +  02J{Y} 

=     rM*^I"1[(D]TC-Mj{-xg  -  axg}  (3.23) 

where 

aMJ     =     C^w-J  (3.24) 

The  terms  on  the  right  side  of  Eq.   3.23  can  be  written  as 

r-M*^_1[(})]T[-Mj{-xg  -  axg}     =     {PF}(-xg  -  axg)  (3.25) 

where 

{PF>     =     ^M^^M^Mjd}  (3.26) 

and  PF.   is  the  participation  factor  for  the  ith  mode.     Substituting  Eq.   3.26 
into  Eq.   3.23 

{¥}  +  r'2Sso)J{Y}  +  F26ho)J{Y}  +  [>2J{Y} 

=     {PF}(-xg  -  axg)  (3.27) 
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For  the  ith  mode,  then, 

Y.   +23     co..Y.   +  23u  to.-Y.   +  o>?.Y.     =     PF.(-x     -ax)  (3.28) 

i  s.   n   i  h.   n   i         n   i  i       g  g 

11  3  3 

Since  the  left  side  of  Eq.  3.28  is  linear  with  respect  to  the  unknown  Y,  the 
maximum  modal  response  in  the  ith  mode  can  be  obtained  by  determining  the 
maximum  response  of  the  single  degree  of  freedom  equation, 

yi  ♦  2^0,,..^  +  ZB^a,,*,  ♦  ^.y.  =  -xg  -  axg  (3.29) 

and  multiplying  this  response  times  the  participation  factor  for  the  ith  mode, 
PF..  The  nonlinear  equations  of  motion  for  the  multi -degree  of  freedom  system, 
Eq.  3.1,  have  now  been  transformed  into  a  set  of  linear  simultaneous  differ- 
ential equations  (Eq.  3.5)  by  the  use  of  Eqs.  3.18  and  3.22.  The  linear 
simultaneous  differential  equations  are   transformed  into  nm  uncoupled,  single 
degree  of  freedom  modal  equations  of  the  form  represented  by  Eq.  3.28  in  the 
normal  coordinate  system  Y  by  the  use  of  Eq.  3.2.  The  total  response  of  the 
system  is  the  result  of  a  combination  of  the  various  modal  responses  represent- 
ed by  Eq.  3.28. 

3.2.5  High  Frequency  Modes 

The  existence  of  a  high  frequency  tail  in  the  response  spectrum 
for  a  single  degree  of  freedom  oscillator  was  demonstrated  and  discussed  in 
Chapter  2.  Further,  the  linearization  of  the  nonlinear  hydrodynamic  drag 
force  was  shown  to  break  down  when  this  high  frequency  tail  began  to  appear. 
A  method  for  determining  the  parameters  necessary  to  estimate  the  magnitude 
of  this  tail  for  each  mode  in  the  multi -degree  of  freedom  system  is  now 
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developed.  Consider  the  equations  of  motion  for  an  undamped  multi -degree  of 
freedom  system 

r*Mj{Z)  +  [K]{Z>  =  -FMj{xg}  -  rKDJ{xs|xsl>        (3.30) 

If  the  system  is  so  stiff  that  {x  }  s  (x  „},  Eq.  3.30  can  be  written  as 

s  g 

C^MJ{Z}  +  [K]{Z>     =     -^Mj{xg}  -  rKDJ   (xg|xg|}  (3.31) 

Transform  Eq.   3.31   into  the  normal   coordinate  system  by  the  use  of  Eq.   3.2 
and  premultiply  by  [<j>] 

[*]TrMj[<J>]{Y}  +  [<J>]T[K][<j>]{Y} 

=     -WTrM^]{I}xg  -  [(f»]TrKD-J{I}(xg|xg|)  (3.32) 

or, 

r-M*^j{Y}  +  rK*j{Y> 

=  -E*]T[^Mj{I}xg  -  [c|»]TrKD-J{I}(xg|xg|)  (3.33) 

Equation  3.33  represents  nm  single  degree  of  freedom  uncoupled  modal  equations 
in  normal  coordinates  for  the  high  frequency  multi-degree  of  freedom  system. 
Premultiply  Eq.  3.33  by  C^M*^]"1,  then 

{Y}  +  [>2J{Y}  =  -{PF}xg  - '{KJj}(xg|xg|)  (3.34) 


where 


{K*}  =  C-M^'^iDflMCnJII}  (3.35) 


For  the  ith  mode  then, 
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KD. 

Y.  +  coT.Y.  =  PF.  (-x  -  of1  x  |x  I)  (3.36) 

i    n  i      i  v  g    PFi  g1  g" 

Again,  since  the  left  side  of  Eq.  3.36  is  linear  with  respect  to  the  unknown 
Y,  the  maximum  modal  response  in  the  ith  mode  can  be  obtained  by  determining 
the  maximum  response  of 

KD. 

y.  +  o)2.y.  =  -x  -  ■==!  x  |x  I  (3.37) 

Ji         nJi      g    PFi   g1  g1 

and  multiplying  this  response  times  the  participation  factor  for  the  ith  mode. 

It  was  shown  in  Chapter  2  that  the  quantity  ly. |maw  can  be  estimated  for 

i  max 


Eq.  3.37  as 


y  • 
lJ-\  'max 


PF. 

l      /  •     \  2 


2  *      g  max 


\  wiiMii 


(x„  maJ  0.38) 


and 


Y.l     =  (PF.)(|y.|   )  (3.39) 

l'max    v  -r^i'max' 


3.2.6  Superposition  of  Modal  Responses 

The  maximum  response  in  any  given  normal  mode  of  the  multi -degree 
of  freedom  system  can  be  computed  from  Eq.  3.28.  However,  the  maximum  re- 
sponse in  one  mode  does  not  necessarily  occur  at  the  same  time  as  the  maximum 
response  in  another  mode,  and  the  time  of  occurrence  of  the  maximum  response 
in  the  time  history  of  the  modal  response  is  not  normally  known.  The  actual 
combination  of  modal  responses  in  a  multi -degree  of  freedom  system  is  very 
complex  and  few  generalities  can  be  made  about  the  combination. 


69 


An  approximation  to  the  actual  maximum  response  of  a  multi -degree 
of  freedom  structure  can  be  made  by  determining  the  maximum  response  in  each 
mode,  and  then  combining  the  maximum  modal  responses  in  some  fashion.  An 
upper  bound  on  the  maximum  response  of  the  structure  can  be  established  by 
assuming  that  the  maximum  modal  responses  occur  at  the  same  time  in  each 
mode  and  are  additive  in  nature.  In  other  words,  the  summation  of  the  ab- 
solute values  of  the  maximum  response  in  each  mode  yields  an  upper  bound  to 
the  response  of  the  structure.  Another  commonly  used  method  of  combining  the 
maximum  modal  responses  is  to  take  the  square  root  of  the  sum  of  the  squares 
of  the  modal  responses.  Other  methods  of  combining  modal  responses  are  some- 
times used.22  The  two  most  commonly  used  methods  are  the  summation  of  the 
absolute  values  of  the  modal  maxima  (ZABS)  and  the  square  root  of  the  sum 
of  the  squares  of  the  modal  maxima  (SRSS).  These  two  methods  are  the  only 
ones  used  in  this  thesis  for  illustrative  purposes.  Tables  3.28  through 
3.33  show  the  results  of  the  superposition  of  maximum  modal  responses  applied 
to  three  of  the  five  degree  of  freedom  systems.  The  drag  force  is  linearized 
in  accordance  with  Eqs.  3.18  and  3.22  (ground  velocity  linearization),  and 
the  amount  of  linear  hydraulic  drag  $.  in  each  mode  is  computed  from  Eq.  3.24. 
The  maximum  modal  response  in  each  mode  is  computed  from  the  response  of  the 
single  degree  of  freedom  modal  equation  (Eq.  3.29),  and  multiplied  times  the 
participation  factor  for  the  mode.  The  SABS  combination  is  seen  to  be  higher 
than  the  true  response  (as  should  be  the  case)  and  the  SRSS  combination  is 
generally  less  than  the  true  response.  No  generalizations  as  to  the  relative 
accuracy  for  different  distributions  of  mass,  stiffness,  and  drag  factor  can 
be  made,  however.13  Thus,  it  is  seen  that,  if  the  maximum  responses  of  the 


70 


single  degree  of  freedom  modal  equations,  Eq.  3.28,  can  be  computed,  the  maxi- 
mum response  of  the  structure  can  be  approximated  by  a  combination  of  the 
modal  maxima. 

3.3  Idealized  Elastic  Deformation  Response  Spectra 

It  was  seen  in  the  preceding  section  that  the  maximum  response  of 
a  multi -degree  of  freedom  system  can  be  approximated  by  combining  maximum 
modal  responses  in  some  fashion.  If  the  maximum  modal  responses  can  be  esti- 
mated in  a  way  that  accounts  for  the  effects  of  the  nonlinear  hydrodynamic 
drag  forces  in  the  system,  a  rapid  approximation  of  the  response  of  the  sys- 
tem can  be  obtained.  The  principles  developed  in  Chapter  2  for  the  single 
degree  of  freedom  oscillator  and  the  method  of  linearizing  the  drag  forces 
developed  in  this  chapter  can  be  combined  to  construct  an  idealized  elastic 
deformation  response  spectrum  for  each  mode  of  the  multi -degree  of  freedom 
system  subjected  to  strong  motion  earthquake  excitation.  The  procedure  for 
the  construction  of  the  modal  elastic  deformation  response  spectra  is  now  set 
forth. 

1.  Compute  the  natural  frequencies,  ["wj,  and  mode  shapes,  [<|>], 
for  the  system. 

2.  Establish  the  structural  damping  ratio,  $s.,  for  each  mode. 
Values  of  3c.  between  0.005  and  0.05  are  reasonable  to  assume 
for  the  type  of  structure  being  investigated  in  this  thesis. 

3.  Determine  x     ,  x     ,  x     ,  and 

g  max   g  max   g  max 

•2,-  , 

g1  g1 

•2 
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for  the  earthquake  to  be  used.  If  the  spectra  are  to  be  con- 
structed for  a  hypothetical  or  "design"  earthquake,  one-half 

of  the  value  of  x     is  a  reasonable  value  to  assume  for  the 
g  max 

quantity 

g'  g 

<> 

Compute  M..,  Kq  ,  PF.,  a,  3n.,  and  the  high  frequency  tail 
estimate  from  the  following  expressions: 


FmVj     =     W^MJW 


(k*}  =  ^^Wr^m 


{PF}     =     rM*^]-1[({)]T[^Mj{I} 


(3.35) 
(3.26) 


/  i 


a     = 


J-l      D.1.1 


nn 


I       M.. 
\  j=l       JJ 


(3.22) 


r23ho)J     =     a[-IJ 


(3.24) 


High  Frequency  Tail 
Estimate  for  the  ith  = 
Mode 


TJ. 
l 


PF 


l 


2  * 
uf.M. 


|   / 


^xg  max^ 


(3.38) 


5. 


With  the  quantities  to..,  $s.,  3h.»  and  the  high  frequency  tail 
estimate,  construct  the  idealized  elastic  deformation  response 
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spectrum  for  the  ith  mode  exactly  as  described  in  Section  2.5 

for  the  single  degree  of  freedom  system. 
It  is  noted  that  the  idealized  response  spectrum  constructed  as 
outlined  above  accounts  for  the  fluid-structure  interaction  in  a  slightly 
conservative  manner.  When  an  idealized  spectrum  has  been  drawn  for  each  mode 
of  the  system,  the  maximum  modal  responses  are  obtained  from  the  spectra, 
multiplied  times  the  appropriate  modal  participation  factor,  and  combined  in 
some  manner  to  yield  an  approximate  maximum  response  of  the  system. 

3.4  Hydrodynamic  Parameters 

An  approximate  method  for  computing  hydrodynamic  parameters  for  a 
multi-degree  of  freedom  system  by  the  superposition  of  normal  mode  responses 
is  developed  in  this  section.  As  explained  in  Chapter  2,  these  hydrodynamic 
parameters  are  absolute  displacement  and  velocity  responses  of  the  mass  (each 
node  in  the  multi-degree  of  freedom  system),  whose  magnitudes  are  desirable 
to  know  in  order  that  the  empirical  drag  and  inertia  coefficients,  Cq  and  Cj, 
may  be  estimated  in  a  reasonable  way.  The  hydrodynamic  parameter  spectrum 
amplification  factors  developed  in  Chapter  2  may  be  used  to  construct  idealized 
spectra  from  which  the  response  quantity  of  interest  in  each  mode  can  be  esti- 
mated. 

Consider  again  the  linear  equations  of  motion, 

r*Hj{Z}  +  [Cs]{Z>  +  [K]{Z}  =  -FMj{xg}  -  rKDLJ{xs)     (3.4) 

Making  the  substitution  {Z}  =  {x  }  -  {x  },  Eq.  3.4  can  be  written  as 

s     g 
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FMj{xs}  +  [Cs]{xs>  +  rKD|_J{xs>  +  [K]{xs> 

=  [K]{xg>  +  [Cs]{xg>  (3.40) 

If  the  velocity  force,  [C  ]{x  },  can  be  dropped  from  Eq.  3.40  for  values 
of  [C  ]  up  to  some  magnitude  without  unreasonable  loss  of  accuracy  in  the 
solution  (this  is  discussed  in  subsequent  paragraphs  of  the  section),  and 
the  drag  factor  matrix  ["*  IC.  J  is  made  to  be  proportional  to  the  mass  matrix 
[""Mj,  then  Eq.  3.40  can  be  written  as 

FMJ{XS)  +  [Cs]{xs>  +  a[-Mj{x$}  +  [K]{xs>  =  [K]{xg}    (3.41) 

In  order  to  transform  Eq.  3.41  into  a  "normal"  coordinate  system,  the  trans- 
formation 


{xs>  =  [4>]{Q}  (3.42) 

is  made.  Substituting  3.42  into  Eq.  3.41  and  premulti plying  the  result  by 

of, 

=     [4>]T.[K]{x„}  (3.43) 


Define 


g 


(Y>     =     l>]T[K]{I}  (3.44) 


and  Eq.   3.43  can  be  written  as 


Fm\J{Q}  +  rC*J{Q}  +  arM*^J{Q}  +  r~K*J{Q}     =     {Y}xn 
*>  y 


(3.45) 
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Equation  3.45  represents  nm  uncoupled  single  degree  of  freedom  modal  equations 
of  motion.  The  equation  of  motion  for  the  ith  mode  is 

M*.Q.     +     C*     Q.     +     aM*.Q.     +     K*.Q.     =     y.x  (3.46) 

nvi  s..j   i  nxi  n   i  'i  g  v 

or, 

<J1  +  2%«A  +  Z\«iA  +  «fal     ■     (^  ><<4iXg>  (3-47) 

Since  the  left  side  of  Eq.  3.47  is  linear  with  respect  to  the  unknown  Q,  the 
response  quantities  of  interest  in  the  ith  mode  can  be  obtained  by  computing 
the  desired  response  quantities  from 

9  9 

q     +     26  coq     +  23hcoq     +     oo  q     =     oi  x  (3.48) 

and  multiplying  the  result  times  the  ith  mode  participation  factor  for  hydro - 
dynamic  parameters.  The  participation  factor  for  hydrodynamic  parameters 
can  be  seen  from  Eqs.  3.44  and  3.47  to  be 

{HDPF}  =  FK^LclOVin}  (3.49) 

Once  the  modal  quantities  are  obtained  in  the  Q  coordinate  system,  Eq.  3.42 
can  be  used  to  obtain  an  approximation  of  the  actual  hydrodynamic  parameter 
of  interest  at  each  node  in  much  the  same  manner  as  an  estimate  of  the  maxi- 
mum deformation  of  each  node  is  obtained  by  modal  superposition. 

To  answer  the  first  question  raised  in  this  section,  the  suitability 
of  dropping  the  velocity  force,  [C  ]{x  },  from  Eq.  3.40,  response  calculations 
were  made  for  the  following  two  equations  of  motion: 

r-Mj{xs>  +  [Cs]{xs>  +  rKDL-J{x$}  +  [K]{xs> 

=  [K]{xg}  +  [Cs]{xg}  (3.40) 
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and 

^Mj{xs>  +  [Cs]{xs>  +  rKDL^]{xs}  +  [K]{x$}     =     [K]{xg}  (3.50) 

The  quantities  {|xcl        },   {<|x    l>},   {|xl       >,  and  {<|x    |>}  for  systems  1, 

S  mdX         b  b  MldA  b 

4,  and  7  in  the  mid-frequency  range  were  calculated  for  the  El  Centro  earth- 
quake from  both  Eqs.  3.40  and  3.50.  The  damping  matrices  [C  ]  used  in  the 
series  of  calculations  were  such  that  they  yielded  damping  ratios  of  .02, 
.05,  and  .10  in  each  mode  and  ["~KDL-J  was  set  to  zero.  For  each  structural 
system  and  value  of  [C  ],  the  hydrodynamic  parameters  computed  from  Eq.  3.40 
were  within  5  percent  of  the  corresponding  response  quantities  computed  from 
Eq.  3.50.  The  results  of  the  calculations  for  system  4  and  a  damping  matrix 
which  yields  .02  in  each  mode  are  shown  in  Table  3.34.  In  another  set  of 
calculations,  the  damping  matrices  [C  ]  were  made  to  be  proportional  to  the 

mass  matrix  and  the  numerical  values  of  [C  ]  were  such  that  the  first  mode 

s 

damping  ratios  were  .02,  .05,  and  .10.  Again,  for  each  structural  system 
and  value  of  [C  ],  the  response  quantities  computed  from  Eq.  3.40  were  within 
5  percent  of  the  corresponding  response  quantities  computed  from  Eq.  3.50. 
The  numerical  results  were  such  that  no  clear  trend  resulted  that  indicated 
which  equation  yielded  consistently  higher  or  lower  results.  The  hydrodynamic 
parameters  were  computed  for  system  4M  in  the  mid-frequency  range,  with  a 
structural  damping  matrix  which  yielded  modal  damping  ratios  of  .02  in  each 
mode,  for  the  El  Centro  earthquake.  The  drag  factor  matrix  was  linearized  by 
using  the  ground  velocity  characteristics  and  made  to  be  proportional  to  the 
mass  matrix.  In  one  set  of  computations,  the  velocity  force  [C  ]{x  }  in  Eq. 
3.40  was  included  and  in  another  set  of  computations,  the  velocity  force 
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[C  ]{x  }  was  neglected.  The  two  sets  of  numerical  results  differed  by  less 
s   g 

than  5  percent. 

From  these  calculations,  it  is  concluded  that  the  velocity  force 
[C  ]{x  }  can  be  dropped  from  Eq.  3.40  without  sacrificing  the  accuracy  of 
the  results  by  more  than  about  5  percent,  for  values  of  [C  ]  which  yield 
modal  damping  ratios  of  up  to  about  .10. 

In  order  to  show  the  similarity  of  form  between  Eq.  3.48  and  the 
linear  equation  of  motion  for  the  single  degree  of  freedom  oscillator  dis- 
cussed  in  Chapter  2,  it  is  instructive  to  recall  Eq.  2.7,  repeated  here  for 
convenience. 

•      *         o 

Z  +  23  oaZ  +  2B.u)Z  +  a)  Z  =  -x  -  23uo)X  (3.51) 

s      n  g    h  g 

Since  Z  =  x  -  x  ,  Eq.  3.51  can  be  written 
s    g 

SL  +  23  cox     +  23ucdx^  +  cA       =     w2x     +  23  wx  (3.52) 

ssshs  s  gsg 

It  was  seen  in  Chapter  2  that,  for  values  of  3total  up  to  about  .10,  the 
hydrodynamic  parameters  computed  from  Eq.  3.51  were  dependent  upon  the  amount 
of  total  damping  in  the  system  (3  +  3h)  and  not  on  the  distribution  of  the 
damping  between  structural  and  hydraulic.  If  all  of  the  damping  in  Eq.  3.51 
is  considered  to  be  structural  and  all  of  the  damping  in  Eq.  3.52  is  considered 
to  be  hydraulic,  the  velocity  force  on  the  right  side  of  both  equations  can 
be  dropped,  and  Eqs.  3.51  and  3.52  become 

Z  +  23tooZ  +  o)2Z  =  -x  (3.51a) 


and 


2       2 

xs  +  23twxs  +  a)  xs  =  (ox  (3.52a) 
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In  other  words,  for  values  of  Bt  up  to  about  .10,  the  hydrodynamic  parameters 

lx     »  <|x  |>,  |x  I   ,  and  <|x  l>  computed  from  Eq.  3.51a  are  within  2 
's'max   '  s1   '  s'max      '  s1     v  ^ 

to  4  percent  of  those  which  would  be  computed  from  Eq.  3.52a  for  the  same 
system  and  base  disturbance. 

It  is  seen  that  Eq.  3.48  is  of  the  same  form  as  Eq.  3.52a,  so  if 
the  hydrodynamic  parameters  can  be  calculated,  or  estimated  (as  from  an  ideal- 
ized response  spectrum)  for  Eq.  3.52a  with  satisfactory  accuracy,  then  modal 

hydrodynamic  parameters  Iql   ,<|q|>,  Iql   ,  and  <|q|>  can  also  be  calculated 
j         j  r  '^'max   Ml   M,max      '^' 

or  estimated.  Once  the  modal  response  quantities  Mmax»  <hl>»  l^lmax»  anc' 
<|q|>  are  obtained  and  multiplied  times  the  appropriate  modal  hydrodynamic 
parameter  participation  factor  to  obtain  |Q|   ,  <|Q|>,  |Q|   ,  and  <|Q|>, 

II  Id  A  IMdA 

the  problem  arises  of  how  to  combine  the  modal  responses. 

To  obtain  an  estimate  of  (|x  I   >  and  {|x  I   ),  the  modal  responses 

s  max        s  max 

can  be  combined  in  any  one  of  several  methods  as  is  the  case  for  obtaining 

{|ZL,„K  An  upper  bound  is  obtained  by  summing  the  absolute  values  of  the 

max 

modal  responses  (ZABS). 

An  upper  bound  on  the  quantities  {<|x  |>}  and  {<|x  |>}  can  also  be 
obtained  by  the  EABS  method  of  modal  combination,  but  this  is  not  readily 
apparent.  That  the  ZABS  method  of  modal  combination  yields  an  upper  bound 
on  the  quantity  {<|x  |>}  (or  {<|x  |>})  is  shown  as  follows: 

{xs>  =  [<i>]{Q}  (3.42) 

Assume  that  the  vector  {Q(t)>  has  been  obtained,  then,  at  the  jth  node  of 
the  system, 

x  (t)  =  ♦J1Q1(t)  +*j2Q2(t)  ♦  .  .  .  *j;nmQnm(t)        (3.53) 
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Taking  the  absolute  value  of  both  sides  of  Eq.  3.53 


|x.  (t) 
J 


(J^Q^t)  +  <J>.2Q2(t)  + 


cj>.   Q   (t) 


(3.54) 


An  upper  bound  on  |x  (t)|  is  seen  from  Eq.  3.54  to  be 

J 


|x    (t)|    <     |<fril  |  IQ-, (t) |  +  ki2MQ2(t)|  + 


k     iiq   (t) 

1Tj  ,nm'  |Xnmv    ' 


Integrating  both  sides  of  Eq.   3.55  over  the  time  of  response, 


(3.55) 


ft  ,      ft 

|xs  (t)|dt    <    j         Cl*j-,  I IQ-, (t) I     +     k-2MQ2(t) 


>•       I  IQ     (t)|]  dt 
vj ,nm' |Xnmv    '  ' J 


(3.56) 


or, 


<|xs  |>    <     k-pMQ^    +     kj2l<|Q2l> 

J 


+  .    .    .    |<J>.        |<|Q     |> 
Irj,nm'    '  xnm' 


(3.57) 


Which  is  seen  to  be  the  IABS  method  of  modal  combination. 

Rules  for  constructing  idealized  hydrodynamic  parameter  response 
spectra  from  which  |q|   ,  <|q|>,  |q|   »  and  <|q|>  may  be  estimated  for  each 

•  • 

mode  (and  then  |Q|   »  <|Q|>»  iQl   •>   and  <|Q|>  are  computed  from  the  q's  and 
|N|max   |X|   |X|max      '  '        v  ^ 

the  hydrodynamic  parameter  participation  factors)  of  the  multi-degree  of 
freedom  system  are  discussed  in  the  next  section. 
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3.5  Idealized  Elastic  Hydrodynamic  Parameter  Response  Spectra 


Idealized  elastic  hydrodynamic  parameter  response  spectra  can  be 
constructed  for  each  mode  of  the  multi -degree  of  freedom  system  in  the  fol- 
lowing manner: 

1.  Compute  the  natural  frequencies,  ["w-J,  and  mode  shapes,  [<J>], 
for  the  system. 

2.  Establish  the  structural  damping  ratio,  3  »  for  each  mode. 

si 

3.  Compute  3.   and  HDPF.  from 


a     = 


nn  \ 

i    V. 

J  =  1         JJ 
nn 

J-l     MjJ      / 


<x2|x   |>  \ 

g1  g1 

•2 
<x  > 

g 


(3.22) 


r23hw-J     =     aMJ 

{HDPF}     =     rKtj_1[(})]T[K]{I} 


(3.24) 
(3.49) 


4.  With  the  quantities  go..,  3<-.,  3h  »  x 


n 


g  max 


(or  <|x  |>),  x 


g  max 


(or  <|x  |>),  and  x     (or  <|x  |>),  construct  the  idealized 
v    '  g1       g  max  v    '  g 

elastic  hydrodynamic  parameter  response  spectrum  for  the  ith 
mode  exactly  as  described  in  Section  2.7  for  the  single  degree 
of  freedom  system. 
When  an  idealized  spectrum  has  been  drawn  for  each  mode  of  the  sys- 
tem, maximum  modal  hydrodynamic  parameters  are  obtained  from  the  spectra, 
multiplied  times  the  appropriate  modal  hydrodynamic  parameter  participation 
factor,  and  combined  in  some  manner  to  yield  an  approximate  solution  of  the 
hydrodynamic  parameter  of  interest. 
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CHAPTER  4 

APPLICATIONS 

4.1  General 

In  Chapters  2  and  3  a  method  was  developed  for  linearizing  the  non- 
linear hydrodynamic  drag  force  on  a  single  degree  of  freedom  oscillator  and 
a  multi-degree  of  freedom  lumped  parameter  model.  The  method  has  been  shown 
to  yield  satisfactory  results  for  systems  subjected  to  earthquake  excitation 
of  the  type  used  in  this  thesis.  Although  the  method  of  linearizing  employed 
in  the  multi -degree  of  freedom  system  places  a  linear  damper  at  the  level  of 
the  top  mass  where  a  nonlinear  damper  does  not  exist,  it  has  been  shown  to 
yield  maximum  response  quantities  which  are  within  a  wery   few  percent  of  the 
response  quantities  computed  for  the  nonlinear  system.  Further,  linearizing 
the  drag  force  by  using  velocity  characteristics  of  the  ground  motion  has  been 
shown  to  yield  generally  conservative  results  in  all  but  the  low  frequency 
range.  The  linearization  technique  has  been  developed  in  such  a  manner  that 
the  resulting  linear  simultaneous  differential  equations  of  motion  for  the 
multi-degree  of  freedom  system  possess  classical  normal  modes  of  response. 
The  existence  of  classical  normal  modes  permits  the  use  of  normal  mode 
superposition,  with  the  modal  response  quantities  determined  by  response 
spectrum  techniques,  as  an  approximate  method  of  computing  the  response  of 
multi -degree  of  freedom  systems.  The  approximate  method  can  be  used  for  dis- 
placement response  quantities  and  hydrodynamic  parameter  response  quantities. 
The  research  reported  in  this  thesis  is  not  intended  to  reinvent  the  method  of 
solution  by  the  superposition  of  normal  mode  solutions  as  determined  by 


81 


response  spectra,  but  is  devoted  primarily  to  developing  the  linearizing 
technique  which  insures  the  existence  of  classical  normal  modes  and  yields 
generally  conservative  results. 

This  chapter  demonstrates  the  use  of  the  principles  and  methods 
developed  and  discussed  in  Chapters  2  and  3  in  computing  the  response  of 
a  series  of  actual  offshore  tower  designs  and  one  of  the  test  systems  des- 
cribed in  Chapter  3. 

4.2  Displacement  Response 

Table  4.1  lists  the  properties  of  a  series  of  offshore  towers  de- 
signed for  water  depths  of  400,  600,  800,  and  1000  feet.1*  A  structural  damp- 
ing matrix  [C  ]  which  yields  modal  damping  ratios  of  .01  in  each  mode  is  com- 
puted for  each  tower  as  described  in  Chapter  3.  The  structures  are  subjected 
to  the  El  Centro  earthquake  depicted  in  Fig.  1.2,  and  various  response  quan- 
tities are  computed.  In  the  first  set  of  computations,  the  equations  of  mo- 
tion are  numerically  integrated  and  response  quantities  are  found  for  the 
case  in  which  there  are  no  hydrodynamic  drag  forces  acting  on  the  structure. 
The  same  response  quantities  are  next  computed  for  the  case  in  which  the  actual 
nonlinear  drag  forces  are  acting  on  the  structure.  Finally,  the  nonlinear 
drag  forces  are  linearized  by  using  the  velocity  characteristics  of  the 
ground  motion,  and  the  response  quantities  are  computed  by  the  combination 
(ZABS  and  SRSS)  of  maximum  modal  responses  found  from  the  true  response  spec- 
trum for  each  of  the  first  five  modes  of  the  structures.  The  amount  of  linear 
hydraulic  damping  in  each  mode  found  from  the  linearization  procedure  is  shown 
in  Table  4.2.  The  response  quantities  of  interest  in  these  calculations  are 
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the  maximum  displacement  between  the  top  deck  and  the  ground,  the  maximum  over- 
turning moment,  and  the  maximum  base  shear  expressed  as  a  percent  of  the  total 
tower  weight,  including  the  added  mass  of  the  structure.  The  results  of  the 
calculations  are  shown  in  Fig.  4.1. 

A  study  of  Fig.  4.1  reveals  that,  for  this  series  of  towers  sub- 
jected to  the  El  Centro  earthquake,  the  nonlinear  hydrodynamic  drag  forces 
reduce  the  response  of  the  tallest  tower  in  the  top  deck  displacement  and 
overturning  moment  quantities  by  approximately  30  percent  from  the  "no  drag" 
case.  The  maximum  base  shear  for  the  tallest  tower  is  reduced  only  about  8 
percent  by  the  nonlinear  hydrodynamic  drag  forces.  Exact  response  calculations 
were  made  for  the  other  three  earthquake  records  used  in  this  thesis  and  the 
reduction  in  response  due  to  drag  was  found  to  be  somewhat  lower  in  magnitude 
than  for  the  El  Centro  earthquake.  For  the  tower  in  1000  feet  of  water,  the 
reduction  (due  to  hydrodynamic  drag)  in  the  maximum  top  deck  displacement 
relative  to  the  ground  was  found  to  be  about  8.5  percent  for  the  Taft  earth- 
quake, 11.5  percent  for  the  Ferndale  earthquake,  and  17  percent  for  the  San 
Fernando  earthquake.  In  general,  the  taller,  more  flexible  towers  showed 
more  reduction  in  response  due  to  the  hydraulic  drag  than  did  the  shorter, 
stiffer  towers.  No  generalizations  can  be  made,  however,  regarding  the  mag- 
nitude of  the  reduction  in  response  relative  to  tower  height  or  stiffness 
since  the  actual  response  is  such  a  complex  combination  of  variables  in  both 
the  structure  and  the  earthquake. 

In  the  superposition  of  maximum  modal  responses  computed  from  the 
true  response  spectrum  for  each  of  the  first  five  modes  of  vibration,  it  was 
found  that  the  higher  modes  contributed  rather  strongly  to  the  maximum  base 
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shear.  The  contribution  of  each  mode  to  the  maximum  base  shear  computed  by 
the  ZABS  combination  of  modal  responses  is  shown  in  Table  4.3.  Although  the 
SRSS  combination  of  modal  responses  appears  from  Fig.  4.1  to  provide  a  very 
good  estimate  of  the  response  of  this  series  of  towers  to  the  El  Centro  earth- 
quake, generalities  concerning  the  accuracy  of  one  method  of  modal  combination 
relative  to  another  method  of  modal  combination  cannot  be  made.13 

The  response  of  the  towers  is  computed  by  the  combination  of  maxi- 
mum modal  responses  (ZABS  and  SRSS)  found  from  the  idealized  deformation  re- 
sponse spectrum  for  each  of  the  first  five  modes  of  vibration.  The  idealized 
spectra  were  constructed  in  accordance  with  the  procedure  described  in  Chapter 
3  for  ground  motion  quantities  equal  to  those  of  the  El  Centro  earthquake.  The 
idealized  spectra  for  the  first  five  modes  of  tower  C  are  shown  in  Fig.  4.2. 
The  estimated  hydraulic  damping  ratios  in  the  3rd,  4th,  and  5th  modes  of  the 
structure  are  such  that  the  resulting  idealized  spectra  for  those  modes  cannot 
be  shown  separately  and  distinctly  on  the  scale  of  Fig.  4.2.  It  is  noted  that, 
for  the  idealized  spectra  shown  in  Fig.  4.2,  the  frequencies  of  the  system  are 
well  removed  from  the  high  frequency  tail  region  of  their  respective  response 
spectrum.  This  was  found  to  be  the  case  for  each  of  the  towers  A,  B,  C,  and 
D.  The  results  of  the  computations  using  the  idealized  spectra  are  shown  in 
Fig.  4.3. 

4.3  Hydrodynamic  Parameters 

To  demonstrate  the  computation  of  the  hydrodynamic  parameters 
lxeLav>  <|x  l>,  |x  I  .  .  and  <|x  l>  at  each  node  of  a  multi-degree  of  freedom 

o  HldX       5        S  ITlaX  S 

system,  the  five  degree  of  freedom  system  4M  is  used.  The  properties  of  this 
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system,  with  a  structural  damping  matrix  that  yields  a  damping  ratio  of  .02 
in  each  mode,  are  shown  in  Table  4.4.  The  magnitude  of  [K]  was  chosen  to 
yield  frequencies  of  the  system  that  fall  into  the  velocity  amplification 
region  of  the  spectrum.  The  hydrodynamic  parameter  response  quantities  are 
computed  by  the  numerical  integration  of  the  equations  of  motion  and  by  the 
combination  (ZABS  and  SRSS)  of  maximum  modal  responses  found  from  the  true 
spectral  response  in  each  mode.  The  results  of  the  computations  are  shown 
in  Table  4.5.  It  is  seen  from  Table  4.5  that  the  combination  of  modal  re- 
sponses provides  results  which  are,  in  general,  sufficiently  close  to  the 
exact  response  at  each  node  to  enable  the  selection  of  inertia  and  drag 
coefficients  C,  and  Cn  with  reasonable  confidence.  In  lieu  of  using  the 
true  spectral  response  quantities,  which  require  the  use  of  a  digital  com- 
puter to  obtain  for  each  mode,  an  idealized  spectrum  can  be  constructed  for 
each  mode  in  the  manner  set  forth  in  Chapter  3. 
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CHAPTER  5 
CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 
5.1  Conclusions 

In  this  thesis,  an  approximate  method  has  been  developed  for  com- 
puting the  response  of  a  multi-degree  of  freedom  system,  partially  submerged 
in  water,  to  earthquake  excitation.  In  Chapter  2,  the  effects  of  the  non- 
linear hydrodynamic  drag  force  on  a  single  degree  of  freedom  oscillator  are 
shown,  and  a  generally  conservative  method  of  linearizing  the  drag  force  is 
developed.  Amplification  factors  are  developed  which  permit  the  construction 
of  idealized  response  spectra  for  deformation  and  hydrodynamic  parameters. 
Chapter  3  extends  the  principles  from  the  single  degree  of  freedom  oscillator 
to  multi-degree  of  freedom  systems  and  develops  a  generally  conservative 
method  of  linearizing  the  drag  forces.  The  method  of  linearizing  insures  the 
existence  of  classical  normal  modes  in  the  equations  of  motion,  which  in  turn 
permits  the  response  quantities  of  interest  to  be  obtained  by  the  superposition 
of  normal  mode  responses.  In  Chapter  4  the  response  to  the  El  Centro  earthquake 
of  a  series  of  offshore  towers  is  computed  using  the  methods  developed  in 
Chapter  3. 

Several  conclusions  may  be  drawn  from  the  development  of  principles 
and  methods  and  the  results  of  sample  problems  in  the  preceding  three  chapters. 

1.  The  maximum  relative  displacement  response  of  a  single  degree 

of  freedom  oscillator,  with  drag  characteristics  described  by 

K  x 
the  parameter  9  max  ,  approaches  a  constant  value  at  high 

frequencies  as  a  result  of  the  nonlinear  hydrodynamic  drag 
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force  growing  very   large  relative  to  the  inertia  force  in  the 
equation  of  motion.  The  magnitude  of  this  relatively  constant 
high  frequency  tail  can  be  estimated  for  earthquake  motions  by 
the  expression 


'  K,  x      '  *2 


'max 


D  g  max 
K 


/  x- 
g  max 

xg  max 


2.  Linearizing  the  drag  force  on  a  single  degree  of  freedom  oscil- 
lator by  using  the  absolute  velocity  response  of  the  oscillator 
yields  slightly  conservative  results  for  the  maximum  response 
to  earthquake  motions  until  the  high  frequency,  high  drag  tail 
begins. 

3.  Linearizing  the  drag  force  on  a  single  degree  of  freedom  oscil- 
lator by  using  the  velocity  characteristics  of  the  ground  motion 
yields  conservative  results  for  the  maximum  response  to  earth- 
quake motions  except  for  low  frequency  systems  and  very  high 
frequency  systems. 

4.  The  conservatism  introduced  by  using  the  ground  motion  velocity 
characteristics  in  the  linearization  of  the  drag  force  is  re- 
duced as  the  amount  of  structural  damping  in  the  system  is  in- 
creased. 

5.  The  maximum  response  of  the  single  degree  of  freedom  oscillator 
to  earthquake  motions  is  sensitive  to  the  total  amount  of  linear, 
velocity  dependent  damping  in  the  system,  and  not  to  the  distri- 
bution of  this  damping  between  structural  and  hydraulic,  up  to 
linear  hydraulic  damping  ratios  of  10  percent  of  critical. 
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6.  The  nonlinear  hydrodynamic  drag  has  \/ery   little  effect  on  the 
elasto-plastic  deformation  response  spectrum  in  the  low  fre- 
quency region. 

7.  A  high  frequency  tail  exists  in  the  elasto-plastic  deformation 
response  spectrum  for  a  system  with  nonlinear  hydrodynamic  drag 
forces.  Large  amounts  of  drag  cause  greater  reductions  in  the 
yield  deformation,  Z  ,  from  the  maximum  elastic  deformation  for 
a  given  ductility  factor  than  do  small  amounts  of  drag. 

8.  The  nonlinear  hydrodynamic  drag  force  has  a  damping  effect  on 
the  elasto-plastic  deformation  response  spectrum  in  the  mid- 
frequency  range.  The  damping  effect  of  the  drag  becomes  less 
pronounced  as  the  ductility  factor  increases. 

9.  Idealized  response  spectra  which  consider  the  effects  of  fluid- 
structure  interaction  can  be  constructed  as  described  in  this 
thesis. 

10.  The  linearization  of  the  drag  forces  on  a  multi -degree  of  freedom 
system  by  using  the  absolute  velocity  response  at  each  node  and 
forcing  the  linear  drag  factor  matrix  [^IC..^]  to  be  proportional 
to  the  mass  matrix  ["Mj  as  described  herein  is  satisfactory. 

11.  Linearizing  the  drag  forces  on  a  multi-degree  of  freedom  system 
by  using  the  velocity  characteristics  of  the  ground  motion  and 
forcing  the  drag  factor  matrix  to  be  proportional  to  the  mass 
matrix  as  described  herein  yields  generally  conservative  results 
except  for  low  frequency  systems. 

12.  For  the  multi-degree  of  freedom  systems  studied  in  this  thesis, 
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the  response  to  earthquake  excitation  is  relatively  insensitive 
to  the  spatial  distribution  of  the  linear  drag  forces. 

13.  The  hydrodynamic  drag  forces  can  cause  a  significant  reduction 
in  the  response  of  lightly  damped  systems  to  earthquake  motions 

14.  The  amount  of  reduction  in  response  due  to  hydrodynamic  drag 
varies  considerably  from  one  response  quantity  to  another  and 
from  one  earthquake  to  another. 

15.  The  reduction  in  response  due  to  hydrodynamic  drag  is  decreased 
as  the  amount  of  structural  damping  in  the  system  is  increased. 

16.  Higher  modes  of  a  multi-degree  of  freedom  system  can  contribute 
appreciably  to  the  total  response  of  the  system. 

5.2  Recommendations  for  Further  Study 

The  approximate  methods  of  analysis  developed  in  this  thesis  are 
limited  by  the  types  of  multi -degree  of  freedom  systems  used  in  the  investi- 
gation and  by  the  inherent  limitations  in  the  hydrodynamic  force  equation 
used.  Further,  the  study  in  this  thesis  is  devoted  to  excitation  of  the 
structure  in  still  water,  and  the  ocean  is  never  completely  at  rest.  There 
are  several  areas  in  which  further  studies  should  be  made  of  the  response  of 
offshore  structures  to  earthquake  excitation. 

1.  The  validity  of  the  lumped  parameter  model  should  be  tested 
through  the  use  of  field  measurements  of  the  response  of  actual 
offshore  structures  to  strong  motion  earthquakes. 

2.  Model  tests  of  single  degree  of  freedom  oscillators  and  multi- 
degree  of  freedom  structures  should  be  conducted  for  the  types 
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of  base  motion  considered  in  this  thesis.  Such  model  tests  will 
be  useful  in  verifying  the  principles  developed  herein. 

3.  Since  the  ocean  is  never  actually  at  rest,  analytical  and  ex- 
perimental studies  should  be  conducted  on  the  response  of  sub- 
merged and  partially  submerged  structures  subjected  to  earthquake 
motions  combined  with  forces  resulting  from  waves  and  currents. 

4.  The  studies  in  this  thesis  assume  ground  motion  in  one  direction 
only.  It  is  recommended  that  studies  be  made  of  the  response 
when  the  two  orthogonal  components  of  the  horizontal  ground  motion 
are  acting  on  the  structure  at  the  same  time. 

5.  The  earthquake  response  of  the  structure  in  torsion  and  in  the 
vertical  direction  should  be  investigated. 

6.  It  is  recommended  that  the  effects  of  foundation  and  soil  con- 
ditions be  included  in  the  model.  There  are  several  reasonable 
approaches  to  modeling  the  foundation  which  should  be  investi- 
gated.23"25 

7.  The  hydrodynamic  force  equation  used  in  this  thesis  (Morison's 
equation)  is  based  on  studies  with  relatively  rigid  piles.  There 
is  some  evidence19'20  that  the  forces  on  flexible,  oscillating 
structures  may  be  better  expressed  through  the  use  of  different 
drag  and  inertia  coefficients,  or  even  through  the  use  of  a 
different  force  expression.  Much  more  information  is  needed 
concerning  the  hydrodynamic  forces  acting  on  flexible,  oscillating 
structures. 

8.  The  linearizing  technique  should  be  tested  for  suitability  of 
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use  when  the  structure  has  radically  different  distribution  of 
mass,  stiffness,  and  drag  properties  from  those  test  structures 
used  in  this  thesis. 
9.  It  is  recommended  that  the  effects  of  the  hydrodynamic  forces 
on  the  inelastic  response  of  the  structure  be  studied  in  greater 
depth.  While  design  for  elastic  response  is  satisfactory  for 
earthquakes  which  may  be  expected  to  occur  frequently,  elastic 
design  for  infrequent,  severe  earthquakes  may  result  in  ex- 
cessively conservative  designs.  Some  response  into  the  in- 
elastic range  should  be  permitted  for  the  occurrence  of  an  un- 
usually strong  earthquake,  and  approximate  methods  are  desirable 
for  the  estimation  of  this  inelastic  response. 
10.  It  is  recommended  that  approximate  methods  be  developed  for 
the  estimation  of  the  response  to  stochastic  representations 
of  the  forces  resulting  from  earthquakes,  waves,  and  currents. 
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Table  2.1 


Characteristics  of  Series  Representation 
of  Base  Acceleration 


Quantity    Maximum  Amplitude    Time  of  Occurrence 
x  2.56  g  0.35  t. 

g  d 

x  3.85  ft/sec  0.42  t, 

9  d 

x  0.20  ft  0.50  t, 

9  d 

x     =  -0.20401   cos  2-rrt  -   6.87108  cos  4iTt  +  38.35016  cos  6TTt 

g 

-  41.35223  cos  8-rrt  -   2.84232  cos  lOirt 
+  12.91937  cos   12-irt 

(x     in  ft/sec2) 
y 
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Table  2.3 

Distribution  of  Total  Damping  Between  Structural 
and  Hydraulic  for  Comparison  Purposes 


Bh  For 

Btotal 

3s  =  0 

3s 

=  .02 

3S  =  .05 

3  =  .10 

s 

.02 

.02 

0 

.05 

.05 

0 

.07 

.05 

.02 

.10 

.10 

.05 

0 

.12 

.10 

.02 

.15 

.10 

.05 

Table  2.4 
Elastic  Deformation  Spectrum  Amplification  Factors 


Ampl- 

ification 

Factor 

etotal^ 

Displacement 

Velocity 

Acceleration 

0 

2.5 

4.0 

6.4 

0.5 

2.2 

3.6 

5.8 

1 

2.0 

3.2 

5.2 

2 

1.8 

2.8 

4.3 

5 

1.4 

1.9 

2.6 

7 

1.2 

1.5 

1.9 

10 

1.1 

1.3 

1.5 

20 

1.0 

1.1 

1.2 
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Table  2.5 

Elastic  Hydrodynamic  Parameter  Spectrum 
Amplification  Factors 


IxJ         Spectrum 
1   s'max     K 

Amplification  Factors 

total  Displacement        Velocity 

(%) 


0  2.4  5.2 

1  2.2  4.4 

2  2.1  3.6 
5  1.8  3.1 

10  1.6  2.1 


<|x  |>  Spectrum 

Amplification  Factors 

total  Displacement        Velocity 

(%) 


0  3.8  10.9 

1  3.4  8.5 

2  3.0  6.8 
5  2.3  4.2 

10  1.8  3.3 
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Table  2.5  (Continued) 


•  • 

1  s 
Ampli 

Lax  sPect™m 
fication  Factors 

3total 
(%) 

Displacement 

Veloci 

ty 

Acceleration 

0 

2.2 

4.1 

7.1 

1 

1.9 

3.1 

4.9 

2 

1.6 

2.6 

3.7 

5 

1.4 

2.0 

3.1 

10 

1.2 

1.8 

2.6 

<|x  |>  Spectrum 
Amplification  Factors 
total    Displacement    Velocity    Acceleration 


0 

3.4 

10.1 

13.5 

1 

3.0 

6.3 

6.2 

2 

2.7 

5.1 

4.3 

5 

2.1 

3.0 

2.9 

10 

1.6 

2.2 

2.2 
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Table  2.6 


Actual  and  Approximate  Reduction  Factors  for  the 
High  Frequency  Region  of  Inelastic  Response 
Spectra.  El  Centro  Earthquake 
K„x 


Ps  "  u.. 

K 

-  U.  JU 

Actual 

Ductility 

y 

Factor 

Z 

y 

(inches) 

'  'max 
(inches) 

Reduction  Factor 
forZy 

y 

2y  -  1 

1 

0.63676 

0.63676 

1.00 

1.00 

1.5 

0.49526 

0.74289 

0.77 

0.75 

2 

0.42835 

0.85670 

0.67 

0.67 

3 

0.35975 

1.07925 

0.56 

0.60 

4 

0.31523 

1.26093 

0.49 

0.57 

Note:  All  calculations  are  at  a  frequency  of  10  cps. 
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Hydrodynamic  Parameters  for  System  4,  Mid-Frequency  Range, 


El   Centro 

Earthquake. 

3     =   .02  in 
s 

Each  Mode 

Parameter 

Node 

Including 

[csKxg} 

Excluding 
CCs]{xg} 

*lxsLax* 

1 

18.95250 

18.98743 

(inches) 

2 

16.31281 

16.37181 

3 

13.45125 

13.47894 

4 

12.30791 

12.21299 

5 

10.53108 

10.50990 

«|xsl» 

1 

6.28175 

6.26993 

(inches) 

2 

5.21131 

5.20904 

3 

4.48454 

4.48830 

4 

3.75600 

3.75798 

5 

2.81634 

2.81841 

{l*sU} 

1 

32.53683 

32.50833 

(in. /sec) 

2 

23.84686 

24.27509 

3 

19.37517 

19.68547 

4 

21.22523 

21.28960 

5 

15.38559 

15.47496 

{<|xs|» 

1 

10.26027 

10.26663 

(in. /sec) 

2 

7.95916 

7.94259 

3 

6.83778 

6.82866 

4 

5.61321 

5.60034 

5 

4.52049 

4.51665 
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Table  4.2 

Estimated  Linear  Hydraulic  Damping  Ratios 
for  Towers  A,  B,  C,  and  D 


Towe 

rs 

Mode 

A 

B 

C 

D 

1 

0.00587 

0.00632 

0.01178 

0.01015 

2 

0.00251 

0.00304 

0.00576 

0.00541 

3 

0.00145 

0.00168 

0.00315 

0.00328 

4 

0.00107 

0.00126 

0.00246 

0.00233 

5 

0.00085 

0.00115 

0.00199 

0.00183 

Table  4.3 

Contributions  of  Each  Mode  to  the  Maximum 
Base  Shear.  ZABS  Modal  Combination 


Contri 

bution 

(as 

%  of 

Total )  of  Mode 

Tower 

1 

2 

3 

4                  5 

A 

32.8 

40.2 

18.7 

6.8             1.5 

B 

33.0 

28.3 

24.0 

13.6             1.1 

C 

19.3 

24.6 

15.6 

35.0             5.5 

D 

38.9 

41.8 

17.9 

1.4             0.0 

Table  4.4 

Properties  of  System  4M,  Mid-Frequency  Range. 
3     =   .02  in  Each  Mode 
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[K]: 


[Cs]: 


.00457 


3  lb-sec  /in 


00457         .00914         .01371       lb-sec2/in.2 


15.32 

-15.32 

0 

0 

0 

-15.32 

45.96 

-30.64 

0 

0 

61.29 

-30.64 

0 

lb/in. 

Sym. 

61.29 

-30.64 
76.61 

0.13917 

-0.06907 

-0.01734 

-0.01149 

-0.00479 

0.24264 

-0.09762 

-0.02638 

-0.00830 

0.28693 

-0.10783 

-0.01633 

lb-sec/in 

Sym. 

0.40007 

-0.12933 
0.58443 

Modal   Properties 


de 

Frequency 
(cps) 

en 

(estimated) 
(none) 

Hydrodynamic  Parameter 
Participation  Factor 
(none) 

1 

0.24705 

0.00794 

2.49115 

2 

0.55706 

0.00352 

1.13695 

3 

0.84337 

0.00233 

0.66289 

4 

1.06231 

0.00185 

0.25179 

5 

1.52207 

0.00129 

0.01573 

Table  4.4  (Continued) 
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Spectral  Response  Quantities  (True  Spectra) 


Mode 


'^'max 
(in.) 

12.618 
13.382 
9.316 
10.376 
10.811 


<|q|> 

(in.) 

4.044 
3.197 
3.096 
2.926 
2.712 


16.648 
22.807 
25.950 
38.248 
31.994 


<|q|> 


1 n 'max 
(in. /sec)    (in. /sec) 


5.295 
7.056 
10.182 
9.798 
9.387 
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Table  4.5 

Hydrodynamic  Parameters  for  System  4M. 
El  Centro  Earthquake 


Response 
Quantity 

Node 

Exact 
Value 

Modal 
Using  True 

ZABS 

Superposition  - 
1  Spectral  Response 

SRSS 

1 

17.82343 

31.32392 

21.00747 

2 

16.05878 

21.29961 

16.04786 

{|x  |    } 
s  max 

3 

13.26795 

18.78246 

13.53510 

(inches) 

4 

11.71001 

17.13259 

11.14007 

5 

10.23945 

11.85761 

6.68999 

1 

5.09378 

9.28465 

6.43887 

2 

4.08737 

6.63971 

5.12691 

(<|x  |>} 
s 

3 

3.44525 

5.81354 

4.30800 

(inches) 

4 

2.63410 

4.96109 

3.32816 

5 

1.45048 

3.42885 

1.91468 

1 

31.62296 

50.63429 

30.31177 

2 

22.94491 

35.58812 

22.48838 

C|x  1   } 

s  max 

3 

17.93618 

30.87149 

19.23872 

(in. /sec) 

4 

19.32336 

26.62459 

16.66792 

5 

13.95569 

22.63737 

12.06787 

1 

8.98920 

16.32071 

9.69329 

2 

6.82424 

11.41066 

7.21779 

{<|xsl» 

3 

5.89449 

10.20436 

6.32494 

(in. /sec) 

4 

4.96801 

8.89230 

5.18860 

5 

4.10293 

7.49280 

4.08744 
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Fig.  1.3  Ground  Acceleration,  Velocity,  and  Displacement  of 
El  Centro  Earthquake,  1940,  N-S  Component 
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Fig.  1.4  Ground  Acceleration,  Velocity,  and  Displacement  of 
Taft  Earthquake,  1952,  N-21-E  Component 
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Fig.  1.5  Ground  Acceleration,  Velocity,  and  Displacement  of 
Ferndale  Earthquake,  1954,  N-44-E  Component 
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Fig.  1.6  Ground  Acceleration,  Velocity,  and  Displacement  of  San 
Fernando  Earthquake  (Pacoima  Dam  Record),  1971, 
S-16-E  Component 
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Fig.  2.21  Elastic  Deformation  Response  Spectra  for  El  Centro  Earthquake 
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(°)  6S  =  5%     6h  =  0,  1,  2,  5,  10* 
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Fig.    2.21    (Continued) 
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Fig.  2.22  Elastic  Deformation  Response  Spectra  for  Taft  Earthquake 
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Fig.    2.22  (Continued) 
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Fig.  2.23  Elastic  Deformation  Response  Spectra  for  Ferndale  Earthquake 
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Fig.  2.23   (Continued) 
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Fig.  2.24  Elastic  Deformation  Response  Spectra  for  San  Fernando  Earthquake 
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Fig.  2.24  (Continued) 
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